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PREFACE 


This is a condensed version of Part I of a two part final report on 
T ASA Grant NGR 19-001-016, Evaluation of the Energy Transfer in the Char 
Zone during Ablation. This part describes the theoretical and experimental 
results obtained for the energy absorbed by the nonequilibrium, equilibrium 
and frozen flow of pyrolysis gases in the char zone for heat shield surface 
temperatures up to 3000°F. In Part II of this report results will be 
presented for the analysis of the energy absorbed in the char zone and the 
decomposition zone during ablation for heat shield surface temperatures up 
to 6000°F. 

Part I of the final report also serves as the Ph.D. dissertation of 
Gary C. April. Reference to the unabridged version of the report is suggested 
for detailed discussion on all phases of the research conducted. 
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ABSTRACT 


The energy transfer associated with the reacting flow of pyrolysis products 
through the char layer of a low density nylon- phenolic resin charring ablator 
was studied experimentally and theoretically. It was found that a non- equilib- 
rium flow model, employing finite reaction rate data for the important reactions 
among the pyrolysis products, was necessary to accurately describe the energy 
transport within the char. The important reactions and kinetic data for a 
temperature range of 500°(533°K) to 3000°F (1925°K), with experimental sim- 
ulation to 2300°F (1535°K), were determined and incorporated into the mathe- 
matical model, called the TEMPRE System, This model, in conjunction with 
experimental results obtained in a Char Zone Thermal Environment Simulator, 
were used to clearly show the shortcomings of the limiting cases of frozen 
and equilibrium flow in predicting the true behavior within the qhar layer. 

A comparison of the experimental data obtained using low density, nylon- 
phenolic resin chars was made with the results obtained using graphite as a 
simulated char. The non-equilibrium flow model accurately predicted energy 
transport in the graphite medium using the same important reactions and 
kinetic data developed for flow through chars. This conclusion was needed 
to justify the use of graphite for the measurement of carbon deposition 
from methane and phenol, and, in the catalyst evaluation studies. 

Carbon deposition and decomposition product distributions were deter- 
mined for methane and phenol using carbon- 14 tracers. The product distri- 
butions were helpful in providing additional evidence that the chemical 
reactions included in the model were correct. The identified products of 
methane and phenol thermal decomposition were carbon monoxide, carbon dioxide, 
methane, ethylene, acetylene and phepol. Carbon deposition measurements 
within the char layer were used to locate the temperature where chemical 
reactions among the pyrolysis products became significant. In general, 
deposition was greatest near the front surface where the temperature varied 
between 1800 - 2300°F (1252 - 1535°K). 

The above results were also used in the catalyst evaluation studies. 

The introduction of a catalyst into the pyrolysis product stream (homogeneous) 
or as a coating on the graphite (heterogeneous) was made to accelerate chemical 
reactions by making them occur at lower temperatures. This resulted in 
higher energy absorption by the pyrolysis products. Bromine (homogeneous) 
had an excellent activity for this by lowering the temperature at which 
reactions start from 1900°F (1312°K) to 1500°F (1084°K) . This effect was 
also measured by comparing the experimental exit gas compositions with the 
compositions calculated by the non- equilibrium flow model, using conventional 
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non-catalytic, kinetic data. The relative difference in these values was a 
good measure of the extent of chemical reaction resulting from the addition of 
the catalyst into the system. 

Molybdenum and tungsten co-catalysts (heterogeneous) had essentially no 
effect in accelerating chemical reactions within the char. The slightly 
different carbon deposition profiles obtained were not within experimental 
accuracy to conclusively indicate a beneficial increase in the rate of 
chemical reactions. Platinum catalyst was tested in earlier experiments with 
no success in accelerating the chemical reactions. 

In addition, the air oxidation of nylon- phenolic resin chars was studied 
to determine the rate of oxidation of the char with distance from the front 
surface. The maximum rate was obtained with air flowing from the heated 
front surface through the char and leaving the rear surface. There was no 
flow of pyrolysis products. With a front surface temperature of 2047°F 
(1395°K) and an air mass flux of 0.035 lb/ft2-sec (0.15 kg/m^-sec), an 81% 
conversion of oxygen was obtained. The gas leaving the back surface contained 
4.0% oxygen. This indicated that oxidation was taking place at all depths 
within the char, and this was confirmed with a non- equilibrium flow calculation. 
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INTRODUCTION 


The Nature of Aerodynamic Heating During Planetary Reentry 

One of the most serious problems encountered when space vehicles reenter 
a planetary atmosphere is aerodynamic heating. Typical reentry velocities 
^or various Earth orbital missions are listed in Table 1. Before a manned 
spacecraft can land safely, these speeds must be reduced to conventional air- 
craft speeds. This can be accomplished by applying a reverse thrust or by 
taking advantage of the frictional resistance of the atmosphere. Since the 
return velocity is of the same order of magnitude as the launch velocity the 
reverse thrust method requires the same quantity of fuel for the reentry phase. 
This doubles the fuel requirement for the mission and makes the added weight 
to the system prohibitive. Hence it is more efficient to use the aerodynamic 
braking method to reduce the vehicle speed to a safe level (1,2,3). 


Table 1. Initial Reentry Velocity and Kinetic Energy for Space Vehicles at 

Various Altitudes Above the Earth (1,2) 


Orbital Altitude 
Nautical Miles 

Velocity 

Ft/ Sec (m/Sec) 

Kinetic 
BTU/ Lb 

Energy 

J/,kg 

300 

26,000 

(7910) 

13,500 

(31.5) 

1000 

27,000 

(8210) 

14,600 

(34.1) 

20000 

33,800 

(10350) 

22,800 

(53.2) 

Circumlunar 

36,000 

(11000) 

26,000 

(60.6) 

Venus or Mars 

37,200 

(11390) 

27,600 

(64.2) 


A high speed orbiting vehicle possesses a large amount of kinetic energy 
(K.E. =1/2 mV^) . In aerodynamic braking, this energy is converted to heat as 
the body descends through the resisting atmosphere. In Table 1 the kinetic 
energy (per unit weight) possessed by a vehicle at various orbital altitudes 
above the Earth is also listed. For example, a 5000 pound (2270kg) vehicle 
having an initial reentry velocity of 26,000 feet per second (7910 m/sec) must 
convert 67,500,000 BTU (71,500 J) of kinetic energy to heat. Consider a 
vehicle that is constructed of structural steel having a specific heat of 
0.117 BTUper pound per degree Fahrenheit (4.2 x 10"^ J/kg-°K). If a temperature 
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increase of 1150°P (900°K) is attained, the maximum amount of heat that is 
absorbed is 675,000 BTU (715 J) , or one percent of the total heat generated. 
It is evident from this example that only a small fraction of the thermal 
energy can be permitted to reach the vehicle without causing destructive 
effects. The remaining large fraction (ninety-nine percent) must be trans- 
ferred to the surroundings by the proper selection of the vehicle shape and 
materials of construction (1,3). 

There are two general classifications of body configuration used in 
ballistic reentry design (3, 4, 5, 6): the slender body and blunt body con» 

figurations. The slender body shape causes a minimum aerodynamic drag 
condition with only slight disturbance of the air flow. This shape produces 
a weak, attached shock wave with a large percentage of the heat generated 
being absorbed by the body. The use of slender body configuration is best 
suited to ballistic missile and supersonic flow applications where heating 
loads are experienced for short periods of time. 

Since reentry of space vehicles requires a maximum amount of energy 
transfer to the atmosphere, the blunt or high drag configuration is more 
applicable. A major portion of the energy is absorbed by the air flowing 
between the strong detached shock wave and the vehicle surface and is 
carried away in the wake behind the craft. The shock layer becomes pro- 
gressively hotter during the course of reentry causing dissociation and 
ionization of the air. This results in heat transfer by conduction, con- 
vection and radiation to the surface of the vehicle. Although the heat 
absorbed is a small portion of the heat generated, the relative amount is 
sufficient to produce surface temperatures in excess of 6000°F (3590°K). 
Therefore, a thermal protection system must be employed to protect the 
vehicle from these high temperatures and heat fluxes. 

One possible solution is to provide enough structural mass to safely 
absorb the heat (3,7). However since most metals are poor heat sinks, this 
method would result in extreme weight penalties. 

Transpiration cooling is a second method (8). This technique protects 
the vehicle by injecting a fluid through openings at the body surface into 
the boundary layer. The injected fluid blocks heat transfer into the 
material and maintains a safe temperature at the space cabin wall. This 
method likewise requires additional weight such as equipment to pump and 
regulate the flow of coolant; and as a result, it is prohibitive in manned 
reentry applications. 

Other methods of heat protection exist, i.e. convective, film and 
radiation cooling, but the most successful technique has been ablative 
cooling (7,9). 
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Ablative Thermal Protection of Planetary Reentry Vehicles 

Ablative cooling is similar to heat sink and transpiration cooling in 
method, but it is drastically different in the mechanisms used to achieve the 
desired results. Ablation sacrifices structural stability to preserve thermal 
“©sistivity by melting, vaporizing and/or subliming relatively thin layers of 
the material at the surface. Although absorption of heat by phase change is 
the distinguishing feature of the process, energy dissipation by radiation, 
conduction, convection, transpiration and chemical reaction is likewise 
achieved (7 ,,10). 

Ideally, an ablative material must possess a low thermal conductivity., 
high heat, capacity and large heat of degradation to effectively restrict 
the extreme temperatures to the surface of the vehicle. Success has been 
achieved employing compositions of nylon, phenolic resin, silicon elastomers 
and epoxy resins (5,9,14). 

There are two kinds of ablative protection systems. One is non-charring, 
and the other is charring. A non-charring ablator is one in which the material 
vaporizes into gases and enters the boundary layer counter to the heat flow. 
This counter flow of mass effectively blocks heat transfer into the. material 
and protects the vehicle. Teflon (polytetrafluoroethylene) is one such non- 
charring compound which undergoes chain shortening steps to form small 
polymer units that eventually vaporize. Extensive research with Teflon 
and other non-charring ablative meterials have been reported over a wide 
range, of applications. 

The charring ablator, on the otherhand, has proven to be one of the 
most successful heat shields for reentry heat protection. Being a combination 
of plastics that decompose to a char of porous carbon and low molecular 
weight gases, it protects the vehicle by conduction, convection, plastic 
decomposition, transpiration, endothermic chemical reactions of the pyrolysis 
gases, reradiation from the char surface, and thickening of the boundary 
layer. The charring ablator is conveniently divided into three separate 
zones which include the plastic decomposition zone, the char zone or layer, 
and the boundary layer as shown in Figure 1, 

In the decomposition zone the virgin plastic degrades to char and low 
molecular weight gases. These gases flow through the char zone and undergo 
chemical reactions such as cracking, free radical formation and ionization. 

A very large quantity of heat is absorbed by these predominantly endothermic 
reactions as the gas temperature increases from the decomposition zone to the 
char surface. These hot gases are then injected into the boundary layer 
with additional absorption of heat due to expansion and further chemical 
reactions. 
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Figure 1. Schematic Diagram of the Various Zones Developed 
During Reentry of a Capsule Protected by an 
Ablative Heat Shield (Char Forming Type) . 




Energy Transfer in the Char Zone of a Charring Ablator 

Each of the above regions has been the subject of a sizeable research 
effort, and various types of matematical models to describe the charring 
ablator process have been developed. These include the transient ablation 
analyses of K.ratsch,et _al (16), Kendall, et_al (17) and Swann, et al (18) and 
in more related areas, transpiration cooling studies by Koh and del Casal 
(19,20,21), flow of methane through porous graphite by Clark (22), and 
carbon deposition studies by Weger, et a_l (23,24). In this research 
a better description of the phenomena taking place in the char zone is 
presented. An accurate description is needed of the energy transfer in 
the char layer and the species compositions and fluxes entering the 
boundary layer. At present these variables are evaluated by considering 
the flow to be either frozen (no reaction) or in chemical equilibrium. 
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For frozen flow the lower limit on the energy transfer is computed 
since the energy absorbed by the pyrolysis products is just the change in 
sensible heat as these gases flow through the porous char. This is the 
simplest case to evaluate mathematically, and the classical transpiration 
cooling solution to the energy equation is applicable (25). 

For equilibrium flow in the char zone the upper limit on the energy 
transfer is obtained since chemical reaction rates are infinitely fast, and 
the composition of the pyrolysis products only vary with the temperature as 
predicted by thermodynamic equilibrium. This approximation gives the maximum 
amount of heat that can be absorbed since the reactions occurring are 
predominantly endothermic. The mathematical description of this case is 
more detailed than for frozen flow since an additional term for heat ab- 
sorption by chemical reaction must be included in the energy equation. 

Many investigators feel this model would more accurately describe, the 
actual behavior in the char zone since the reaction rates should be very 
fast at the high temperatures encountered. 

For a more accurate description of the reacting flow in the char zone 
the kinetics of the chemical reactions must be included in solving the 
energy equation. The solution is more complex than the limiting cases because 
compositions of the pyrolysis products must be calculated from the reaction 
rate expressions which are differential equations. Of all the possible 
reactions that could occur in the char zone within the temperature range 
encountered, the ones that actually occur must be selected and included in the 
analysis . 

In addition, experiments must be conducted to assure the theoretical 
model accurately predicts the energy transfer in the char zone. This can 
be accomplished by flowing a mixture of compounds typical of the actual 
pyrolysis gases through chars formed in arc- jet heaters. The chars can be 
radiantly heated to simulate the surface heating during reentry. Gases 
entering and leaving the char zone can be analyzed to determine the extent 
of the reactions taking place in the char. Thus the accuracy of the 
mathematical computations can be assessed. 

Furthermore, the results of the analysis, referred to as the non- 
equilibrium model, can be compared with the limiting cases. In this way 
the limitations incurred by assuming equilibrium or frozen flow are evaluated. 
A detailed investigation into the types of reactions occurring, the amount 
of carbon deposition taking place in the char layer, and ways to make the 
ablation of char forming materials more efficient is determined. 
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MATHEMATICAL ANALYSIS FOR REACTING FLOW IN THE 
CHAR ZONE OF A CHARRING ABLATOR 


General Description 

The momentum, energy and mass transfer associated with the. flow of pyrolysis 
products through the char layer of a char- forming ablative plastic is considered. 
The pyrolysis products, formed by the thermal degradation of the plastic heat 
shield, enter the char layer at the decomposition temperature of the plastic. 

The products experience a temperature increase as they flow through the char 
and undergo thermal cracking to lower molecular weight species which react 
with each other and with the carbonaceous char layer. These predominantly 
endothermic reactions are important modes of energy absorption and must be 
included in any realistic analysis of the energy transfer in the char layer. 

The mathematical model describing the transport phenomena taking 
place has the form of a one dimensional and steady flow. A schematic 
diagram showing the pyrolysis gas flowing through the porous char layer is 
depicted in Figure 2. As indicated, the pyrolysis products enter the char 
at the decomposition temperature, T 0 , and exit at a higher front surface 
temperature, T^. Changes in the mass flux of the various species within 
the char occur as a result of chemical reactions at finite reaction rates, 

R-i . A pressure gradient (P 0 -Pl) across the char is also experienced. 

Restrictions and Assumptions 

The particular restrictions and assumptions used to simplify the equations 
of change are listed in Table 2. Justification of these statements is presented 
in detail in the original work (26) . 



Figure 2. Schematic Diagram of the Char Layer. 
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Table 2. Restrictions and Assumptions Used to Simplify the Equations 
of Change for Reacting Gas Flow in the Char Zone 



Restriction ■- Assumption 

Justification 

i. 

One dimensional flow normal to 
the front surface 

Radius of curuature of capsules 
large compared with thickness 

2. 

Steady flow 

Char thickness remains constant after 
short, initial transient period 

3. 

Ideal pyrolysis gas mixtures 

High temperature, low pressure 

4. 

Variable gas properties 

Large temperature range (500-3000 o F) 

5. 

Constant char properties 
(except char conductivity) 

Carbon deposition small therefore, 
porosity and permeability constant. 
Char thermal conductivity varies 
with temperature. 

6. 

Thermal equilibrium between 
char and pyrolysis gases 

Generally small, especially in 
ablative applications where the 
char back surface temperature is 
nearly equal to the decomposition 
temperature. 

7. 

Momentum transfer by modified 
Darcy's Law 

Inertial effects important 
because of high mass flux values 

8. 

PV work negligible 

Small in comparison to convective and 
conductive transport 

9. 

Negligible viscous dissapation 

Velocity and viscosity of the 
gas mixture are small 

10. 

Biffusional transport omitted 

Sm^ll in comparison with 
bulk fluid transport. 


11. Work against gravity omitted 


Non existent for horizontal flow 



Equations of Change 


The application of the above simplifications to the general equations of 
change result in the important differential equations for describing the flow 
of pyrolysis products in the char layer. These, shown in Table 3, are the 
■pecies continuity equation, the momentum equation, the energy equation and 
the surface heat flux equation. 


Table 3. Summary of the Important Equations Related to the Flow of 
Pyrolysis Products in the Char Zone. 

Species Continuity Eauation: 

|g( Pi v) = Ri ( 1 ) 

Momentum Eauation (Darcy's Law): 

P = tp2 L + 2E[ J L 0) (W) (j*-) dz 

+ J L P(T )(w ) 2 dz"j 1/2 
° % J ( 2 ) 

Energy Equation: 

K 

e 5 p W p H = *z£ k e!i3 "® cRc + i ^ 1 HiR d < 3 ) 

Heat Flux Equation: 

K 

’cz - <lL' q o ) - lSlf L ‘VpiV 1 

i O 

K+l L 

+ S (* H^R^dz ( 4 ) 

i=l "o 


A detailed derivation of these equations from the general equations is 
presented in reference (26) . 










Boundary Conditions 


There are two important sets of boundary conditions that can be used to 
solve the equations of change for flow of reacting gases in the char zone. 

These are shown in Table 4. The first set specifies the pressure and temp- 
rature at the front surface, and, the temperature and pyrolysis gas composition 
entering the back surface of the char. These conditions for mass flux, W, 
as a parameter make the solution of the energy equation a two point boundary 
value problem 1 . 


Table 4. Boundary Conditions for the Solution of the 
Char Zone Equation of Change. 



Iterative Solution 

Non-iterative Solution 

Location 

Boundary Conditions 

Boundary Conditions 

Back 

T = T 0 

T = T 0 

Surface 

(z=o) 

x i = x io 

q p = - (ke -^) 0 

Front 

Surface 

(z=L) 

T = T l 
P = P L 

P = P L 

L (6) 


This requires an iterative solution. Because PV work is negligible in the 
energy equation, the momentum equation (Modified Darcy ' s Law) can be solved for 
the pressure distribution after a solution of the energy equation is obtained 
using an average pressure in the char. The heat flux at the char front surface 
is also calculated using equation (4). The computer program used for these 
solutions is called the Iterative TEMPRE System (IT). 

The second set of boundary conditions specifies the temperature, initial 
pyrolysis gas composition and the sum of the heat of pyrolysis and the heat 
conducted in the virgin plastic at the back surface of the char, q^. For 
mass flux, W, and qp as parameters, the solution of the energy equation as 
an initial value problem can be obtained. This does not require an iterative 
solution as did the first case. To be useful, however, a parametric study of 
W andqp is needed over the range of values expected during reentry. 
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In order to calculate the pressure distribution within the char, the 
front surface pressure is again specified. The heat flux and pressure 
distribution calculations are the same as those made in the Iterative 
TEMPRE System. The computer program for this method of solution is called 
the Non-Iterative TEMPRE System (NIT) . 

Both systems are discussed in detail in reference (26) where a complete 
block flow diagram and program listing are presented. The particular ap- 
plication of the equations of change and the above boundary conditions 
to frozen, equilibrium and nonequilibrium flow of pyrolysis gases within 
the char layer will be developed in the. following sections. 

Application of the Transport Equations to Frozen, Equilibrium 
and Non-Equilibrium Flow in the Char 

As mentioned in the Introduction, there are two limiting cases currently 
used to simplify the analysis of the flow of pyrolysis gases through the 
char zone. These are to consider the flow to be either frozen or in 
thermodynamic equilibrium. This research deals with the development of a 
third model, one for non- equilibrium flow, which will predict more 
accurately the actual behavior within the char layer. In this section 
the equations of continuity, momentum, energy and surface heat transfer 
will be applied to develop each of the three flow models. In a subsequent 
section, the solutions of the particular equations for each model will 
be compared with each other and with experimental data. In this way the 
the accuracy of the analysis for non-equilibrium flow in the char zone 
can be evaluated and the extent to which the two limiting cases predict 
the behavior can be determined. 

Frozen flow. -The frozen flow model is an idealization in which the 
pyrolysis products entering the char zone do not change in composition as they 
pass through the medium. Therefore, any benefit from the energy absorbed by the 
predominently endothermic chemical reactions which occur between the gases 
and char is not obtained. It specifies a lower limit on the amount of energy 
absorbed in the char zone. Of the three cases it is the simplest because 
the chemical reaction terms in the energy and heat flux equations are zero: 

K+l 

S % Ri = 0 (7) 

i=l 

Applying this to the equations of change developed previously results 
in the following simplifications for frozen flow in the char layer. 


Continuity Equation 


Momentum Equation 
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Heat Flux Equation 


( 10 ) 


K T l 

q = (q T - q ) = £ T eW G . X.dT 
4 cz V4 L 4 o' p pi i 

i=I »T_ 


Energy Equation 


SG P W P 


dT = 4 - 

T[z dz 


[ k e 41 ] 


( 11 ) 


The numerical solution of these equations will be discussed in a later 
section. 

Equilibrium flow. -The equilibrium flow in the char zone gives an upper 
limit on the amount of heat that can be absorbed within the char zone. The 
reason is that the reactions occuring within the char are predominately 
endothermic. The set of equations used to describe flow for this case is 
the same as the equations developed previously: continuity (8), momentum 

(9), energy (3), and heat flux (4). The distinguishing feature lies in 
the method used to calculate the energy absorption by chemical reactions: 

K+l 

£ HjRj 0 (12) 

i=l 

Rewriting the species continuity equation, (1), in terms of the mole flux 
of species i gives: 


R i - ■S (pv) ' £ (13> 

K+i 

Therefore, in order to evaluate the term, tke mass flux, W, and 

the mass fraction, x^,‘ of the species in the gas and solid phases must be 
known as a function of temperature. The species composition and molal ratio 
of gases to carbon are a function of temperature, pressure and elemental 
composition of the virgin plastic and can be calculated by one of the many 
approaches in the literature (27,28). In this study, the free energy 
minimization technique was used (26). 

Non-equilibrium flow. -The equations of change for non- equilibrium flow 
are the same as those developed for equilibrium flow. The evaluation of the 
energy absorption by chemical reaction, however, is dependent on the finite 
reaction rates of each reaction occurring within the char zone. Therefore, 
in addition to temperature and pressure, the reaction rate, R^, is also a 
function of the mass flux and composition entering the char layer. This 
requires a knowledge of the specific reactions taking place within the char 
and the associated kinetic data; i.e., frequency factor and activation energy. 
The following paragraphs will present a method for determining the important 
reactions and the technique for using the kinetic data to predict the actual 
behavior within the char zone. 
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Reaction rate and rate constant equation. -In general, a chemical reaction 
can be written in the following form: 


q q 

2 r^ jA^ ^ E pijA^; j = l,2,3,...,m (14) 

i=l i=l 


For this j-th chemical reaction, r^j and are the stoichiometric coefficients 
of the reactants and products respectively ror species A^. There are a total 
of q chemical species and m chemical reactions in the system. 

The rate of reaction of the i-th species, R^, is given by the following 
equation for the m chemical reactions: 


R. 


m 

E (p 
j-1 


ij 


- r , 




)[k 


f h= 


rr c. 


ij-l 


k=l 


rj 


q P ' 

TT C . 

k=l X 


ij 


]; 


i = 1,2,3, ...,q (15) 

where c^ is the concentration of component i and r';[j and p'fj represent 
the power on the concentrations. It is not necessary for these to be equal 
to r-jj and Pi j • The forward and reverse reaction rate constants are k.^ and 
k r j. This is the equation for R^ that is used in the computer implemented 
numerical solution of the transport equations. Equation (15) is a very con- 
venient and general formulation for the reaction rate of the j-th species 
in m simultaneous chemical reactions. The stoichiometric coefficients and 
the powers on the concentrations are each conveniently represented as a 
matrix. 


In addition to the above, the reaction rate constants are also required. 
A general form for the rate constant of the i-th chemical reaction is: 


kj = AjT JExp (-E j/RT) ; j = l,2,3,...,m (16) 

where As is the frequency factor and E., the energy of activation. 

J J 


With these mathematical generalizations the analysis is resolved to one 
of selecting the important chemical reactions taking place in the char zone 
along with precise values of the rate constants. These equations and associated 
rate constants can then be used in equation (15) to calculate the reaction 
rate, R^, which is needed to solve the transport equations. 


The non- equilibrium flow model is an order of magnitude more complex 
than the chemical equilibrium computations. The latter case involves the 
solution of a set of algebraic equations with the energy and momentum equations. 
The former requires the solution of a set of non-linear, ordinary differential 
equations (species continuity equations) simultaneously with the energy 
equation. There is also the additional difficulty of determining all of the 
important chemical reactions that take place in the system and the initial 
composition of pyrolysis products entering the char zone. For equilibrium 
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flow, only the elemental composition is required. Finally, there is the 
laborious task of collecting and evaluating reaction kinetic data appearing 
in the literature. 

Chemical reactions in the char layer. -The general nature of the 
reactions occurring in the char zone of a charring ablator has been qualitatively 
established (29). Between 500°F (533°K) and 3000°F (1925°K), the primary 
reactions occurring are hydrocarbon cracking reactions of high molecular weight 
species to lower molecular weight species (ultimately H 2 , CO 2 , CO, H 2 O, etc.). 

From 3000°F (1925 °K) to about 6000°F (3590°K), free radical and recombination 
reactions take place, with ionization reactions beginning to appear at the 
upper end of this temperature region. This study is concerned with the 
reactions occurring below 3000°F (1925°K); i.e., hydrocarbon cracking reactions, 
primarily. Free radical, recombination and ionization reactions are not 
considered as they do not take place. Fortunately, this greatly simplifies 
the already complex system of reactions that could occur. Establishing the 
fact that reacting flow in the char zone can be accurately described in the 
temperature range from 500° (533°K) to 3000°F (1925°K) establishes the 
basis for extending the investigation to temperatures in the 3000°F (1925°K) to 
6000°F (3590°K) range. 

Pike (30), in an effort to condense the large assortment of reaction 
kinetic data available in the literature, has compiled a detailed listing of 
reactions and the corresponding kinetic data for the C-H-O-N system. In a 
subsequent report (31) this and other kinetic data were analysed. The calculation 
of the isothermal conversion for each reaction possible in the char zone over 
a wide range of temperatures was included. This formed one method for deter- 
mining the important reactions taking place in the char zone as discussed below. 

Criterion for reaction selection. -The rate of a chemical reaction in- 
creases with temperature. For a particular reaction, if a significant conver- 
sion of reactants to products is obtained with the char at a uniform and 
specified temperature, then it can be assumed that there may be a significant 
conversion when a temperature gradient exists in the char layer with the front 
surface at this specified temperature. Thus, this reaction is considered 
important. An example of this behavior is the reaction illustrated in Figure 3 
as a plot of conversion versus temperature. For another reaction if there is 
no conversion in the char at a uniform higher temperature (1925°K), there 
will be no conversion when there is a temperature gradient in the criar. 

Therefore, this reaction can be omitted 'when there is a temperature gradient 
in the char with a front surface temperature of 1925 °K. However, the products 
formed by the reaction may indeed be included as important components subject 
to further reactions. This logic forms the basis of the isothermal analysis 
of reaction kinetics data in the literature (31). 

Although the reactions postulated to take place over a given temperature 
range (533°K to 1925°K) could occur, they are also restricted by the components 
initially present in the pyrolysis gas stream. For example, the reaction of 
butane and oxygen forming carbon dioxide and water is a likely candidate between 
500 F (533 K) and 3000 F (1925°K); however, the absence of either butane or 
oxygen eliminates the reaction from consideration. A considerable amount of 
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caution must be exercised to make certain that the components are not formed 
at some time within the char, making the reaction an important part of the 
isystem. 


The ten chemical reactions 
which accurately describe the 
energy transfer associated with 
reacting pyrolysis gas flow 
through the char zone of a 
charring ablator between 500°F 
(533°K) and 3000°F (1925°K) 
are listed in Table 5. In 
addition to the criterion 
just discussed for eleminating 
reactions, the composition 
of the pyrolysis products 
entering the char at the back 
surface was also considered. 

Composition of pyrolysis 
products. -As discussed above, 
the pyrolysis product composition 
at the back Surface of the char 
is important in reducing the 
total possible reactions 
likely to occur within a 
specified temperature range. 

First attempts to study the 
non-equilibrium flow of 
pyrolysis products through 
a porous char layer relied 
on two separate sources for 
estimating the gas composition 
entering the char at the back 
surface. The first method 
which served as an order of 
magnitude analysis was the 
equilibrium compositions 
calculated by the free energy 
minimization method (26). 

The second was analyses of the degradation products of low density nylon- 
phenolic resin composites by pyrolysis gas chromatography (32, 33). 

The unavailability of accurate analytical procedures and thermo- physical 
data for the high molecular weight pyrolysis products (i.e., phenol, cresol, 
toluene, etc.) left a region of definite uncertainty. As a result, the 
major components of the pyrolysis products were methane, hydrogen, carbon 
dioxide, carbon monoxide and nitrogen by many (34,35) with unknown quantities 
of water and high molecular weight residues completing the analysis. 

Subsequent research by Sykes (36) confirmed the presence of phenol- based 
materials as primary constituents in the high molecular weight residues. 
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Table 5. Important Reactions and Associated Kinetic Data for the Pyrolysis Product 

Species in the Char Zone Between 533 - 1925 K. 


General Form of the Reactions: aA + bB + ... + = rR + aS + ... + 

General Rate Constant Equation: k = k°T“ s Exp(-E/RT) 




Reaction 

Rate Law 

Activation Energy 
E, (KJ/Kg-mole) 

Frequency 

Factor 

CH 4 = 1/2 H 2 + 1/2 C 2 H 

6 k f A 

38.0 

14 

7.6x10 * 

o 2 h 6 - c 2 h 4 +. h 2 

k fA 

28.0 

14 

3. 1x10 * 

C ? H /( = C ? H, + H ? 

kfA 

16.0 

2.6xl0 8 * 

C 2 H 2 = 2C + H 2 

kjA 2 

4.0 

7 

2.1x10 ** 




12 

C + 2H 2 - CH 4 

k f 

7.0 

2.0x10 *** 

C 6 H 6° + H 2 = H 2° + C 6 H 

6 k f A 

18.0 

2,0xl0 13 * 

C 6 H 6 = 3 C 2 H 2 

k^A 

14.0 

1.4xl0 9 * 

C + H 2 0 = CO + H 2 

kfAB 

33.0 

1.2xl0 9 ** 

CO + H 2 0 = H 2 + C0 2 

k f AB 

12.0 

l.OxlcP ** 

C + C0„ = 2 CO 

kfA-k r R 2 

20.0 

l.OxlO 6 * 

z 

24.0 

1.0x10" 12 ** 


s I References 


42,43 

42,50,55 



66.71.60.61, 

52 . 53 . 64.67 

65.60.61.62, 

53 . 64 . 67.68 



* 1 st Order (sec -1 ) ** 2 nd Order (M 3 /Kg-mole-sec) *** 0 th Order (Kg -mole /M 3 ) 















Table 6 represents a more precise pyrolysis gas composition obtained by Sykes (36). 
Very good agreement was obtained by comparing the reported experimental results 
with an overall energy balance technique using heats of formation and heats 
of combustion data for the virgin plastic and the experimentally determined 
pyrolysis product composition. 

The simultated pyrolysis gas composition used in the experimental program is 
also listed in Table 6. 


[ Table 6. Pyrolysis Products 

Resulting From the Thermal Decomposition, of Nylon- 

Phenolic Resin Composites Including the High Molecular Weight Species. 


Sykes (36> Mole % 






(Nylon- Phenolic Resin) 

Sykes (33) 

Simulated 




• ■ Mo l e % 

Pyrolysis 


By Flash 

By Pyrolysis 

(Phenolic 

Product 



Pyrolysis 

at 50°C 

Resin Only) 

Composition 

Component 

at 1073°K 

Increments 







to 900°C 







: fl47.5°K) 





Phenol 

8.0 

7.6 

6.7 


-| 


Methylphenol 

4.0 

4.0 

2.1 




Dimethylphenol 

2.8 

3.2 

0.0 




Trimethylphenol 

2.4 

2.7 

0.0 


6.2 


Benzene 

0.4 

0.4 

0.2 




Toluene 

0.0 

0.2 

0.4 




Cyclopentanone 

. 1.2 

2.5 

0.0 


» — 


Hydrogen 

30.2 

32.2 

47.1 


33.4 


Methane 

3.8 

4.0 

9.4 


6.7 


Carbon Monoxide 

5.3 

3.6 

5.2 


3.7 


Carbon Dioxide 

12.1 

6.8 

1.5 


1.1 


Water 

20.1 

23.3 

22.0 


48.9 


Ammonia 

0.8 

0.4 

0.0 


0.0 


Unidentified 

9.2 

10.1 

5.2 


CLO 


Totals 

100.0 

100.0 

100.0 

100.0 


Physical and thermodynamic properties . -In any real problem where the 
temperature gradient varies over a wide range (>1000°K), changes in the physical 
and thermodynamic properties as a function of temperature occur. For the 
multicomponent flow of a reacting gas within a porous char, composition change 
by chemical reaction is also important. The equations used in this research 
for calculating the variation in physical and thermodynamic properties are 
shown in Table 7. 
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Table 7. Equations Expressing the Physical and Thermodynamic 
Properties as a Function of Temperature 


Gas Phase Properties : 

1. Thermal conductivity: 
of a pure gas component 


of a gas mixture 


2. Viscosity: 

of a pure gas component 


of a gas mixture 


where 


3. Heat capacity 

4. Free energy 


5 . Enthalpy 


Char Properties: 

1. Porosity 

2. Permeability 

3. Thermal Conductivity 


k gi 

= 2.6693 x 10" 5 

m 1/2 

(Cvi+4.47) 

(17) 


<r 2 0v 

k k 


(18) 

kg 

= Ein i kg i / E n£ 
i=l i=l 


^gi 

= 2.6693x10-3 

cr\ 

((Mvi)T) 1 / 2 

(19) 


K K 

, -1-1 

(20) 


= E |J-i j 1+Ei£ 
i-1 L j=l 

JlM- 

5ij(nj/ni)J 


= Q+(p, i /iJ ) p 1 ^ 2 (M w j/M wi ) 1/4 



j^V 8 (l+(M wi /M w j ) 1/2 _1 

(21) 

Cpi/R = 

a t + bj_T + CiT 2 + d t T 3 + ei T 4 

(22) 

O 

& 

H- 

II 

ai (l-log e T) - H T - c iT 2 
^ 6 


RT 

- ii T 3 - fi T 4 - 1± - g . 
12 1 20 1 T i 

(23) 

(H° T )i . 

a i + + fi T 2 + ^t 3 


RT 



+ li T 4 + fi 



5 T 

(24) 

Assumed 

constant and equal to the average 



effective porosity. 

Average values obtained from a Carman - 
Cozeny plot of reported experimental data (26), 

Linear least squares fit of available 
experimental data (26) . 
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Numerical Solution of the Equations of Change 

Prior to the selection of a numerical method, the accuracy of the numerical 
solution desired must be specified. This determines the interval size needed 
in the analysis which effects the round-off errors associated with the calculation 
of the solution. If a relatively small interval size is used, the round-off 
errors may be intolerable. On the other hand, a very large interval size could 
produce large truncation errors which results in a solution that does not approach 
the true solution. One technique used to determine if a specific interval size 
is a reasonable choice involves the computation of the solution for the particular 
interval selected, and interval sizes reduced by a factor of one half of each 
preceeding value. Comparison of the solutions for each interval chosen should 
reveal when the approximate solution approaches the true solution or when 
round-off errors make the calculated solution invalid. 

In general, increased accuracy requires increased complexity in (or order 
of) the numerical method used. Therefore, there is an optimum decision to be 
made between the nearness of the approximate solution to the real solution and 
the computational time required to obtain the solution. The best approach is 
obviously the one that minimizes truncation and round-off errors. 

Numerical solution of the differential energy equation. -The energy equation 
describing the flow of pyrolysis gases through the char zone of a charring ablator 
is a second order, non-linear differential equation with variable coefficients. 

In order to obtain a solution to this equation, a numerical integration technique 
must be used. This requires the equation to be transformed into a finite 
difference form which can be integrated stepwise on a digital computer. Of 
the various methods available, all require the specification of the order and 
stepsize to achieve the desired solution. Two general types of numerical 
integration schemes are commonly used; the self-starting methods, ranging 
from the Euler equation to the Runge-Kutta series, and, the predictor-corrector 
methods which require the specification or calculation of several initial points 
to start the procedure. The number of points is proportional to the order of 
the particular equation used. Because of the accuracy and relatively straight- 
forward nature of the self- starting methods, a fourth order Runge-Kutta formula 
was selected. 

The formulae as they apply to the differential energy equation are presented 
in Table 8 . The truncation error is of the order O(h^) where h is the step size 
(37). The parameters A]_, A 2 , A 3 , and A 4 in equation (28) correspond to the 
Runge-Kutta parameters generated by the numerical integration of the energy 
equation. 

Therefore, in the non- equilibrium flow analysis, the solution of the 
energy and continuity equations are calculated simultaneously because of the 
interdependence of the temperature and the mole flux (composition). 

Numerical solution of the species continuity equation. -In the non- equilibrium 
flow analysis, the species continuity equation expressing the mole (or mass) flux 
of each species as a function of temperature must be solved simultaneously with 
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the differential energy equation. Rewriting equation (13) in terms of the 
mole flux, Nj[, of species i gives: 


dNi R, i 

— = (^-)(— = l —) 
dT M w i dT/dz 


Table 8. Fourth Order Runge-Kutta Formulae for Solving the 
Differential Energy Equation for Flow of Pyrolysis 
Gases Through the Char Zone 


Formulae for the Flow of 


iis Gases Through Chars 


t N+ 1 “ t H + h t T N + t< A 1 + a 2 + A 3^ 


T »1 ■ T i + 7 <A 1 + 2A 2 + ^3 + V 
6 

A x = h(r ) - .£ k e^(T^> + 22 


?] 


" W C„e dk e /dT EH 

A 2 - »(t;+ 1/2 A,) — ^ — (Ti + 1/2 A,) + — 


i^-i 

e _ 


6 J T+-T'+h. 


, . , W a C D e dk e /dT , , EH 

A 3 = h(Tg+l/ 2 A 2 ) -J2-E § (Tg+1/ 2 A 2 ) + — 

k e k e k, 


i^i 

e T n + -T' + -A 
- 1 ° 2 ° 8 2 


a 4 - her; + A 3 ) 


WpCpS 


dk e / dT 


: (t;+a 3 ) + 


c e 


t 0 + hT ; + 2 a 3 


where the temperature gradient, dT/dz, is calculated in the energy equation 
solution. Substitution for R£ from equation (15) results in the final form 
of the equation used to calculate the species mole flux, N is in the TEMPRE system. 
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The solution of the K+l first order differential equations for the mole 
flux, Ni, is obtained by numerical integration using a fourth order Runge- 
Kutta formula (37) shown in Table 9. 


Table 9. Fourth Order Runge-Kutta Formulae for Solving 
The Species Continuity Equation. 



The interdependence of the energy and species continuity equations can be seen 
by noting the presence of the Aj values in the above equations (31). These 
values are the Runge-Kutta parameters calculated for the energy equation 
solution at the same position within the char. A one-to-one correspondence 
between the Aj and j values must exist to obtain a non- equilibrium flow 
solution. To start t&e integration, the initial composition, temperature 
and temperature gradient at the back surface are used to calculate the 
first intermediate temperature and Runge-Kutta parameter, A^ . These 
values are then used to estimate the mole flux of each species at the same 
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intermediate position within the char zone. The intermediate flux values 
are then substituted into the energy equation for the calculation of the 
next intermediate temperature. This procedure is continued until the four 
(for a fourth order analysis) parameters associated with each differential 
equation (continuity and energy) are calculated. At this point the 
temperature and concentration at a position advanced one interval unit 
into the char are calculated using equations (26) and (30) . 

This technique is repeated to the second boundary of the problem. 

If a non-iterative set of boundary conditions is specified (i.e., one 
boundary completely specified), the solution to the energy equation is 
obtained in one cycle and the calculation of the pressure and heat flux 
distributions can be started. If, however, a two point boundary value 
problem is being solved, a guessed value of the gradient must be used to 
start the solution. When the temperature at the front surface is calculated, 
it must be compared with the specified value. A calculated value which is 
too high or low requires an adjustment of the initial gradient and repetition 
of the entire procedure. Therefore, in addition to the Runge-Kutta logic, 
a convergence procedure to approach the known front surface temperature is 
also required. 

Convergence techniques for the iterative Runge-Kutta analysis. - The 
specific convergence scheme used in the Iterative TEMPRE System is sub- 
divided into three main parts. These include: (1) a gross correction 

procedure, (2) a fine correction procedure, and, (3) a limit or check 
procedure. Each method uses a simplified formula which bases the corrected 
adjustment on the relative overshoot or undershoot of the fromt surface 
temperature . 

Gross correction procedure . -This method is used when the absolute 
value of the difference in the calculated and true front surface temp- 
erature exceeds a specified tolerance limit (usually > 57 0 of the specified 
front surface temperature). When this condition is met, adjustment of the 
initial slope is made using the following equation: 


T n 

°new 


t’ (i + x L 

old 


T *• T 

^true ^calculated 


J true 


]) 


(32) 


where X is an adjustable parameter from 0 to 1. 

Fine correction procedure . -This method is employed when the calculated 
front surface temperature falls between tolerance limits of 0.05 > 1 - T/Tl > 
0.01 of the specified value or when two complete cycles have been calculated. 
In this case a more precise adjustment is used taking advantage of the pre- 
vious calculations to rapidly approach the specified front temperature and 
to prevent oscillations around that point. The specific equation used to 
adjust the initial gradient for the next iteration is: 
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T ** T **" 

^true ^calc. this cycle _ 

T T 

_ ^calc. this cycle “ ^calc. last cycle _ 

( 33 ) 

Very good convergence is obtained with these simple equations over a 
range of front surface temperatures between 1500 - 3000°F (1090 - 1925 °K) 
and mass flux values between 0.0003 - 0.1 lb/ft^-sec (0.00015 - 0.5 kg/m~ 
-sec). Three iterations are required to obtain a solution within the stated 
tolerance range. 

Limit or check procedure. -A limit corresponding to a 20% overshoot 
of the front surface temperature is used to terminate the calculations at 
any point within the char where this condition is violated. A reduction 
in the initial gradient is made according to equation (32) and the 
procedure is restarted. 

In addition to the above, several checks are included in the Runge- 
Kutta analysis to insure the calculated and/or adjusted values of the 
gradient remain non-negative and non- zero. The program logic includes 
an adjustment calculation similar to equation (32) which, on repeated 
violations, instructs the system to print pertinent information for 
diagnosis of the problem. With these simple adjustment equations, the 
iterative solution is obtained very quickly and with sufficient accuracy 
to insure a good approximation to the real solution. 

A comparison of the temperature profiles obtained for the frozen 
flow energy equation (11) at various interval sizes is presented in 
Table 10. Based on the presented results, an interval size of one hundred 
units (or 1/400 ths of an inch for a one- quarter thick char) was used. 


t' - t' 

°last cycle ®this cycle 


T 0new T 0 o ld + 


Numerical Solutions of the Heat Flux and Momentum Equations 


The heat flux and momentum equations are solved after a valid temper- 
ature profile is calculated. Both heat flux and pressure are uncoupled 
from the energy equation by considering the energy dissapation by PV 
work small in comparison with other modes of energy transfer in the system. 
The equations for the heat flux and pressure are first order differential 
equations with variable coefficients. 


K T L K+l L 

q = E J eW p Cp x. dT + £ J E^dz 
z i=l T q i=l z 



(4) 

( 2 ) 
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Table 10. Comparison of Various Runge-Kutta Increment Sizes for the Frozen Flow 

Variable Physical Properties Model 


Dimens ionless 


Temperature (°K) 


Char Distance 


Runge-Kutta 

Increment Size 


(z/L) 

10 

50 

100 

400 

0.0 

533.0000 

533.0000 

533.0000 

533.0000 

0.2 

606.0118 

602.7934 

602.7104 

602.6779 

0.4 

713.1568 

710.5278 

710.4730 

710.4393 

0.6 

869.7659 

867.8951 

867.8761 

867.8412 

0.8 

1085.5353 

1084.6029 

1084.5700 

1084.5336 

1.0 

1366.3333 

1366.3333 

1366.3333 

1366.3333 

Conditions : 

W = 0.25 Kg/m^-sec. 

L = 0.0063 M c - 

0.8 

Gas Composition (Mole/Mole Gas)j 




CO » 0.245, 

C0 2 = 0.046, N 2 = 0. 

073, CH 4 = 0.570 

, C.H = 0.068 
o o 






A Simpson's Rule integration technique is used to obtain solutions for 
the integral terms in these equations. The general formula for the 
Simpson's Rule analysis is (37): 


J 2n fdz = i 

2 o 3 


f o + 4 ( f i + f 3 + ••• + f 2n-l) 

+ 2 (f 2 + f 4 + ...+ f 2n . 2 ) + £ 2 1 - £(4) 

J 90 

where f (^) represents the truncation error. 

90 

In terms of the heat flux and pressure equations, the Simpson's 
Rule functions are: 


( 34 ) 


. dz = J- (UO (W) (I-)dz 
P 1 J V ^w 

(35) 

^ dz = J‘^(— ) (W) 2 dz 
P 2 J M„ 

(36) 

4 dz = feW C x.dT 
n l •> P p i 

(37) 

h 2 dz = j‘ H iV z 

(38) 


Once again the step- size is an important parameter that must be 
specified. A similar procedure as that described for the Runge-Kutta 
analysis to obtain the optimum value of the interval size was used. This 
minimizes the errors and maximizes the accuracy of the approximate solution. 
A comparison of the solutions for various step sizes is presented in 
Table 11. An interval size of twenty steps was used in this analysis. 
Calculation of the pressure profile within the char is made from the 
front surface where P = Pp. to the rear surface pressure. Results are 
reported as a pressure distribution and a surface heat flux which 
correspond to the total heat transferred within the char zone. 


Summary of the Theoretical Development of the Equations of Change 

for Flow in the Char Zone 

The equations of change (continuity, momentum, and energy) were 
developed for modeling the reacting flow of pyrolysis gases through the 
char zone of a charring ablator. The particular application to frozen, 
equilibrium, and non- equilibrium flow were discussed along with methods 
for determining the heat absorption by chemical reaction for the latter 
two models. 
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Table H* Comparison o£ Various Simpson's Rule 
Increment Sizes for the Frozen Flow, Variable 
Physical Properties Model 


Dimens ionless 
Char Distance 
(Z/L) 

20 

Pressure 
Simpson's Rule 
50 

( Kg/m 2 ) 
Increment Size 
100 

200 

d.oo 

10877.9605 

10877.9590 

10877.9565 

10877.9560 

0.33 

10867.0735 

10867.0720 

10867.0695 

10867.0690 

0.67 

10843.6865 

10843.6840 

10843.6820 

10843.6820 

1.00 

10800.0000 

10800.0000 

10800.0000 

10800.0000 


Conditions: W = 0.25 Kg/m^-sec L ~ 0.0063 M € 

Gas Composition (Mole/Mole Gas): 

CO = 0.245, C0 2 = 0.046, N 2 = 0.073, CH 4 = 0.570, C^ = 0.068 


i 




Similarly, the momentum and heat flux equations were integrated using 
Simpson’s Rule. The solutions were presented as a pressure distribution 
across the char and a heat flux at the char front surface. A fourth order 
Runge-Kutta integration of the species continuity equation, coupled with 
the energy equation, produced concentration gradients for each pyrolysis 
gas species as a function of the char distance. 

Specific values of the solutions will be presented and compared with 
experimental data obtained on the Char Zone Thermal Environment Simulator 
in a subsequent section. The experimental apparatus and procedure follows. 
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CHAR ZONE THERMAL ENVIRONMENT SIMULATOR 


Experimental Equipment 

The experimental data presented in this study were obtained using an 
apparatus that simulated the flow of pyrolysis gases through the char 
layer of a charring ablator during reentry. Low density nylon-phenolic 
resin chars, graphite and carbon were obtained from the National Aeronautics 
and Space Administration's Entry Structures Branch at the Langley Research 
Center. These specimens were placed in a metal holder with the front 
surface exposed to a bank of infrared quartz lamps used to simulate the 
high temperatures experienced by a reentering vehicle. Simulated pyrolysis 
gases were passed through the char from the rear surface to the heated 
front surface. The exit gases were sampled and analyzed for comparison 
with the inlet gas composition to determine whether chemical reactions had 
pccurred within the char. These results were also compared with the 
Calculated exit gas compositions for frozen, equilibrium and non-equilibrium 
fElow conditions in the TEMPRE system. This was the method used to determine 
the accuracy of the non-equilibrium model for predicting the heat transfer 
for the flow of reacting pyrolysis products in the char zone. A schematic 
diagram of the Char Zone Thermal Environment Simulator is shown in Figure 4. 
A detailed description of the apparatus, including dimensions, materials 
of construction and assembly diagrams are included in reference (26). 
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Experimental Procedure 


The procedure used for all experiments with the Char Zone Thermal 
Environment Simulator considered safety as the prime requisite. The 
recommendations for safe handling of hazardous materials outlined in 
Dangerous Properties of Industrial Materials by N. I. Sax (38) and 
’“Chemical Safety Data Sheets” by the Manufacturing Chemists' Association 
(39, 40, 41) were used as guides. Proper ventilation of the laboratory 
and various safety devices including exhaust hoods, eye wash basin, 
safety shower, fire extinguisher and self-contained breathing equipment 
were available and employed. 

The experimental procedure for the operation of the Char Zone Thermal 
Environment Simulator was divided into five main portions and are summarized 
in Table 12. 


Table 12. 

Experimental Procedure for the Operation of the 

Char Zone Thermal Environment Simulator. j 

j Phase 

Function(s) 

Pre- startup 


Test instrumentation and equipment with helium 
(Argon) flow. 

Startup 


Turn lamps on. Obtain steady state with lamps 
on full power. 

Experimental 

Introduce pyrolysis gases. Record results of | 
steady state runs. I 

Shutdown 


Turn lamps off. Turn pyrolysis gas feed off. 
Introduce helium (Argon) to flush system. 

Analytical 


Analyse gas, liquid and radioactive labeled 
samples . 


The analytical phase was carried out one day after the experiment 
and involved those analyses listed in Table 13. Reference to the original 
research report (26) is made for greater detail and explanation of the 
various procedures summarized in the above paragraphs. 
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Table 13. Typical Analyses of Samples From the Char Zone 
Thermal Environment Simulator 


Gases : H 2 , CO, CH^ Activated Charcoal 

Argon Gas Carrier (30cc./min) 
Temperature 280°F (410°K) 

H 2 , CH^, CO 2 , C 2 H 4 , C 2 H 2 Porapak S 

Argon Carrier (20cc/min) 
Temperature 165°F (350°K) 


Liquids : Phenol in aqueous solutions ASTM Bulletin D 2145 


Radioactive : 

Liquids & Gases : Effluent from gas chromatograph is combusted 

forming CO 2 and H 2 O. CO 2 is absorbed in one 
molar Hyamine Hydroxide solution and analysed 
on a liquid scintillation spectrometer. 

Solids : Chars are combusted in a tubular furnace to 

CO 2 and H 2 O. CO 2 is absorbed in one molar 
Hyamine Hydroxide solution and analysed on 
a liquid scintillation spectrometer. 
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RESULTS OF THE NON- EQUILIBRIUM ANALYSIS OF PYROLYSIS GAS FLOW IN TEE CHAR ZONE 


Comparison of the Non- Equilibrium Flow Results with Experimental 
Data using Low Density Nylon- Phenolic Resin Chars. 

In order to obtain solutions to $he equations of change, the important 
chemical reactions occurring within the char, and the pyrolysis gas 
composition entering at the back surface must be specified. This im- 
formation is summarized in Tables 5 and 6, respectively. The values of the 
activation energy and frequency factor for each reaction is compared to 
literature values in reference (26). In all cases very good agreement was 
obtained. The use of these data in the non-equilibrium flow model allows 
the calculation of the solutions to the energy, momentum and heat flux 
equations describing reacting flow in the char zone. 

The basis for evaluating the non-equilibirum flow model as an accurate 
and realistic analysis of energy transfer in the char layer is by comparison 
of the experimentally determined exit product compositions from the Char j 
Zone Thermal Environment Simulator with the calculated compositions in the 
non-equilibrium flow model. This was done for seventeen experiments using 
low density nylon- phenolic resin chars in which front surface temperatures 
of 1575-2300°F (1130-1535°K) and pyrolysis gas mass flux values of 0.00003- 
0.108 lb/ft^-sec) (0.00015-0.540 kg/m^-sec) were studied. 

In order to obtain a significant conversion of the pyrolysis gases 
over the temperature range investigated, the residence time within the 
porous char was increased by decreasing the mass flux. Thus the combination 
of the highest temperature and lowest mass flux produced a maximum 
conversion resulting from chemical reactions. 

Typical experimental results for the entire range of parameters 
studied are presented in Table 14 . Also listed are the exit gas com- 
positions resulting from the frqzen (or initial composition) equilibrium 
and non- equilibrium calculated results. The excellent agreement of the 
calculated non- equilibrium compositions with the experimentally determined 
exit gas compositions obtained in the Char Zone Thermal Environment Simulator 
prove the analysis to be realistic and accurate over the entire range of 
temperature and mass fluxes investigated. In addition, a definite transition 
from a frozen condition, corresponding to relatively low temperatures and 
high mass flux values, to a non- equilibrium condition, corresponding to 
higher temperature and lower mass fluxes is noted.- Also, in all cases 
except for those in which the mass flux values are extremely small 
(less than 0.00015 kg/m^-sec), the compositions predicted by the 
equilibrium flow model are inaccurate and unrealistic. A more direct 
comparison of the three models (frozen, equilibrium and non- equilibrium) 
will be discussed in the following section with emphasis on the limit- 
ations of the two ideal flow models (frozen and equilibrium). 
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Table 14, Comparison of Calculated and Experimental Exit Gas Compositions from the Char Zone 
Thermal Environment Simulator for Mass Flux Rates Between 0.00655 andi 0.540 Kg/M^-sec 
and Front Surface Temperatures Between 1195 an< j 1545 ®k. 


RUN NUMBER 
MASS FLUX 
FRONT TEMP 
BACK TEMP 

FLOW 

MODEL 

H 2 

MOLE 

% 

CH. 

4 

MOLE 

7. 

CO 

MOLE 

% 

co 2 

MOLE 

% 

n 2 

MOLE 

% 

H 2 ° 

MOLE 

7. 

C,H,0 
6 6 

MOLE 

7. 

C 2«4 

MOLE 

7. 

C 2 H 2 

MOLE 

% 

XVIII-56 

FROZEN 

28.9 

6.4 

3.3 

0.8 

0.0 

53.3 

6.9 

0.0 

Bfi] 

0.01040 

EQUILIBRIUM 

59.7 

3.3 

26.7 

4.0 

0.0 

6.3 

0.0 

0.0 

Kb 

1200°K 

NON- EQUILIBRIUM 

23.0 

11.1 

2.8 

2.1 

0.0 

53.8 

7.1 

0.0 

0.0 

810°K 

EXPERIMENTAL 

24.8 

10.3 

3.4 

1.3 


53.7 

6.5 

0.0 

0.0 


XIX- 60 
0.01240 
1390°K 
1020°K 


XIX- 61 
0.01100 
1450°K 
860°K 





FROZEN 
EQUILIBRIUM 
NON- EQUILIBRIUM 
EXPERIMENTAL 


FROZEN 
EQUILIBRIUM 
NON- EQUILIBRIUM 
EXPERIMENTAL 


FROZEN 
EQUILIBRIUM 
NON- EQUILIBRIUM 
EXPERIMENTAL- 


FROZEN 
EQUILIBRIUM 
NON- EQUILIBRIUM 
EXPERIMENTAL 


FROZEN 
EQUILIBRIUM 
NON- EQUILIBRIUM 
EXPERIMENTAL 


FROZEN 
EQUILIBRIUM 
NON- EQUILIBRIUM 
EXPERIMENTAL 





























































Comparison of the Non-Equilibrium Flow Results with the Results 
of the Frozen and Equilibrium Flow Analyses 

The frozen and equilibrium flow models bracket the non- equilibrium case* 
Frozen flow corresponds to a system in which no chemical reactions occur, 
vjhile equilibrium flow refers to a system of species undergoing chemical 
reactions which are at equilibrium (a function of temperature and pressure 
only). Since the non-equilibrium flow model predicts the actual behavior, 
comparison of the exit gas compositions, temperature and pressure distri- 
butions, and surface heat flux for each model will determine the accuracy 
of the two limited flow analyses in predicting the energy transfer within 
the char layer. 

These results are presented in Figure 5 and Table 15 for a mass flux 
rate of 0.05 lb/ft^- sec (0.25 kg/m^-sec), a front surface temperature of 
1500°F (1090°K)and a back surface temperature of 500°F (533°K) . The char 
porosity is 0.8 and the char thickness is 0.25 inches (0.0064 m) . As 
seen the temperature profile of the non-equilibrium flow analysis is 
identical to the frozen flow temperature profile. The relative closeness 
of these two models is likewise seen by comparing the exit gas composition, 
pressure drop across the char and surface heat flux in Table 15. Therefore, 
at the above conditions there is little evidence of chemical reactions 
in the char and the energy transfer closely predicted by the frozen 
flow model . 

In Figure 6 and Table 16 the same results are presented for a front 
surface temperature of 2000°F (1367°K). Although the reported non-equilibrium 
values are again very nearly equal to the frozen flow results, a noticeable 
change, especially in the concentration profile and surface heat flux is 
observed. This indicates chemical reactions among the species within the 
char layer. 

A continuation of the analysis for a front surface temperature of 
2500°F (1645°K) and 3000°F (1925°K) in Figures 7 and 8 and Tables 17 and 
18, respectively, shows a more dramatic change which is reflected by a down- 
ward shift of the temperature profile toward the equilibrium curve and 
corresponding rapid changes in the concentration profile. Chemical reactions 
are obviously a very important mode of energy contribution under these last 
two sets of conditions. 

It is not possible to extend the analysis to temperatures above 
3000°F (1925 °K) since the chemical behavior within the char will not be 
predicted by the chemical reactions of Table 5 . In this event additional 
reactions must be included to accurately describe the energy transfer 
within the char zone. This extension is discussed by del Valle, et.al. (73). 

In addition to the above comparisons, the temperature profile, 
surface heat flux and pressure drop across the char are compared in 
Figures 7 and 9 for two pyrolysis product compositions. The results 
for the first simulated pyrolysis product composition, which was based 
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0.0 0.25 0.50 0.75 1.0 


Dimensionless Char Distance 

(z/L) 

Figure 5. Temperature Profile for the Frozen, Equilibrium, 
and Non-Equilibrium Flow of Pyrolysis Gases Through the 
Char Zone of a Nylon-Phenolic Resin Ablator. 
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Table 15. Results of the Analyses of Frozen, Equilibriura and Non-Equilibrium Flow of 
Pyrolysis Gas Products Through a One-Quarter Inch Thick Low Density Nylon- 
- Phenolic Resin Ablative Heat Shield Char at 1090°K. 


Conditions: W = 0.25 Kg/M 2 -sec 


T q = 533°K 


Char Distance 

Inlet 

At 

z/L = Q.! 

> 

Flow Model (s) 

(All) 

FF 

EF 

NEF 

Mass Flux (Kg/M2 vo ids“ s ec) 

0.3125 

0.3125 

0.3125 


Temperature (°k) 

533.0 

630.0 

615.0 

630.0 


Composition (Mole% ): 

Hydrogen 

Methane 

Phenol 

Water 

Carbon Monoxide 
Carbon Dioxide 


0.0063 M Feed Mix 4 


At z/L =1.0 


EF 


0.3125 0.3135 0.3125 


33.4 

6.7 

6.2 

48.9 

3.7 


9.2 

33.3 

0.0 

48.4 
0.1 
8.9 


33.4 

6.7 

6.2 

48.9 

3.7 

1.1 


Model 


Pressure Drop Across Char 
Kg/M 2 


Heat Flux at Surface 
KJ/M 2 -sec 


Frozen Flow 
Non-Equilibrium Flow 
Equilibrium Flow 


29.0 

29.5 

26.5 


340.02 

350.35 

509.74 










































an 



(z/L) 


Figure 6. Temperature Profile foT the Frozen, Equilibrium, 
and Non-Equilibrium Flow of Pyrolysis Gases Through the 
Char Zone of a Nylon-Phenolic Resin Ablator. 
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Table 16. Results of the Analyses of Frozen, Equilibrium and Non-Equilibrium Flow of 

Pyrolysis Gas Products Through a One-Quarter Inch Thick Low Density Nylon- 
Phenolic Resin Ablative Heat Shield Char at 1368°K. 

Conditions: W - 0.25 Kg/M 2 - 

sec € 


T o = 

533°K 

L « 0.0063 M 

Feed Mix 4 

Char Distance 

Inlet 

At z/L * 

0.5 

At 

z/L - 1. 

0 

Flow Model (s) 

(All) 

FF 

EF 

NEF 

FF 

EF 

NEF 

Mass Flux (Kg/M 2 volds -sec) 

B9 

0.3125 

0.3135 


m 

0.3140 

0.3140 

Temperature (°{Q 

533.0 

690.4 

650.7 

690.0 

1368.0 

1368.0 

1368.0 

Composition (MoleZ): 
Hydrogen 

33.4 

33.4 

3.2 . 

33.4 

33.4 

18.2 

33.5 

Methane 

6.7 

6.7 

35.3 

6.7 

6.7 

28.8 

6.8 

Phenol 

6.2 

6.2 

0.0 

6.2 

6.2 

0.0 

6.2 

Water 

48.9 

48.9 

54.6 

48.9 

48.9 

42.5 

48.3 

Carbon Monoxide 

3.7 

3.7 

0.0 

3.7 

3.7 

0.3 

4.1 

Carbon Dioxide 

1.1 

1.1 

6.9 

1.1 

1.1 

10.1 

1.2 


Model 


Pressure Drop Across Char 
Kg/M 2 


Heat Flux at Char Surface 
KJ/M 2 -sec 


Frozen Flow 
Non-Equilibrium Flow 
Equilibrium Flow 


36.0 

36.0 

30.0 


550.77 

555.72 

1081.52 





























Temperature ( a F) 



Dimensionless Char Distance 
(s/L) 


Figure 7 . Temperature Profile for the Frozen, Equilibrium, 
and Non-Equilibrium Flow of Pyrolysis Gases Through the 
Char Zone of a Nylon-Phenolic Resin Ablator. 
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Table 17. Results of the Analyses of Frozen, Equilibrium and Non-Equilibrium Flow of 
Pyrolysis Gas Products Through a One-Quarter Inch Thick Low Density Nylon- 
Phenolic Resin Ablative Heat Shield Char at 1650°K. 


Conditions: 


W =0.25 Kg/M Z -sec 


£ = 0 . 8 


T q = 533 K 


= 0.0063 M 


Feed Mix 4 


Char Distance 


Inlet 


At z/L = 0.5 


At z/L =1.0 


Flow Mod el (s) 


(All) 


Mass Flux (Kg/M 2 voids -sec) | 0.3125 0.3125 0.3140 | 0.3125 | 0.3125 0.3070 0.3855 


Temperature (®K) 


Composition (Mole %) : 


Hydrogen 

Methane 

Phenol 

Water 

Carbon Monoxide 
Carbon Dioxide 


Frozen Flow 
Non-Equilibrium Flow 
Equilibrium Flow 


533.0 



690.3 


3.3 

35.4 

0.0 

54.1 

0.0 

7.1 


756.2 I 1650.0 1650.0 


33.4 

6.7 
6.2 

48.9 

3.7 


1650.0 




Pressure Drop Across Char 


Kg/M 2 



Heat Flux, at Char Surface 
KJ/M 2 -sec 


771.65 

1217.48 

2118.71 










































Temperature (®F) 



Dimensionless Char Distance 

<z/L) 


Figure 8. Temperature Profile for the Frozen, Equilibrium, 
and Non-Equilibrium Flow of Pyrolysis Gases Through the 
Char Zone of a Nylon-Phenolic Resin Ablator. 
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Table 18. Results of the Analyses of Frozen, Equilibrium and Non-Equilibr ium Flow of 
Pyrolysis Gas Products Through A One-Quarter Inch Thick Low Density Nylon- 
Phenolic Resin Ablative Heat Shield Char at 1940°K. 

Conditions : W = 0.25 Kg/M 2 -sec £ * 0.8 Tq = 533 °r L = 0.0063 M Feed Mix * 


Char Distance 

Inlet 

At z/L = 0.5 

At z/L =1.0 

Flow Model (s) 

(All) 

FF 

EF 

NEF 

FF 

EF 

NEF 

Mass Flux ( Kg/ni^-yoids " ) 

0.3135 

0.3135 

0.3140 

0.3135 

0.3135 

0.3570 

0.5480 

Temperature (°K) 

533.0 

895.0 

720.3 

785.8 

1940.0 

1940.0 

1940.0 

Composition (Mole %): 

Hydrogen 

Methane 

Phenol 

Water 

Carbon Monoxide 
Carbon Dioxide 

33.4 - 

6.7 
6.2 

48.9 

3.7 
1.1 

33.4 

6.7 
6.2 

48.9 

3.7 

1.1 , 

3.5 

35.4 

0.0 

53.9 

0.0 

7.2 

33.4 

6.7 
6.2 

48.9 

3.7 
1.1 

33.4 

6.7 
6.2 

48.9 

3.7 
JL-1 

63.9 

1.1 

0.0 

2.2 

31.2 

1.6 

o.o (a) 

54.0 

2.4 

0.0 

28.9 

17.4 

Model 

Pressure Drop Across Char 
Kg/M2 

Heat Flux at Char Surface 
KJ/M 2 -sec 

Frozen Flow 
Non-Equilibrium Flow 
Equilibrium Flow 

56.5 
55.0 

36.5 

1017.28 

2126.85 

3690.50 


Note: (a) Non- Equilibrium Flow Model Requires Additional Important Reactions and Associated 


Kinetics Data to Accurately Describe the Energy Transfer Near 1940°C. 
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on experimentally measured and computed equilibrium compositions excluding the 
high molecular weight cyclic compounds, is shown in Figure 9. The results 
in Figure 7 are for the more accurate composition based on pyrolysis gas 
chromatographic analyses of Sykes (33). A mass flux rate of 0.05 Ib/ft^-sec 
(0.25 kg/m^-sec), front and back surface temperature of 2500°F (1644°K) 
and 500°F (533°K), respectively, and a one-quarter inch (0.0064 m) thick low 
density nylon- phenolic resin char (e = 0.8) were the conditions for each 
case presented. 

Comparison of the temperature profiles for frozen and non- equilibrium 
flow show the same overall behavior; i.e., a downward shift by the non- 
equilibrium curves indicating a higher energy absorption due to chemical 
reactions between the pyrolysis products. On the otherhand, a noticeable 
difference in the equilibrium curves is observed. For the more accurate 
pyrolysis gas composition (Figure 7), the characteristically sharp downward 
shift of the equilibrium curve observed in Figure 9 does not occur. 

The explanation for this difference will point out the inadequacy of 
the equilibrium flow model in . predicting the true behavior within the 
char zone. Because the results of Figure 9 were calculated for a pyrolysis 
product composition largely based on equilibrium calculated values, the inlet 
gas composition to the char was very nearly equal to the values calculated 
using the equilibrium flow model. As a result, very little energy absorption 
was omitted from the analysis due to the very small ehange in the inlet 
compositions which were already approximated as an equilibrium composition. 
This resulted in the relatively flat curve over nearly three quarters of 
the total char thickness. 

On the otherhand, the more accurate pyrolysis product composition 
estimated from experimental data and dependent on finite reaction rates 
governing plastic decomposition, is far removed from the equilibrium 
calculated compositions. However, as these concentrations of gases are 
introduced into the equilibrium flow analysis, an abrupt adjustment to the 
equilibrium compositions calculated by minimizing the free energy is 
experienced. This is especially noted for phenol and methane which have 
initial compositions of 6.2 mole % and 6.7 mole %, respectively, but are 
immediately changed to 0.0 mole % and 35.3 mole % at the back surface 
temperature of 500°F (533°K) . This erroneous adjustment results in the 
loss of energy absorption which results from phenol decomposition at 
finite reaction rates and causes the curve to more closely approach the 
frozen and non- equilibrium curves. 

These same conclusions regarding the inability of the equilibrium 
flow model to accurately predict the true behavior is further shown in 
Table 19. Here, the surface heat flux values are compared for each 
model and for each simulated pyrolysis product composition. In comparing 
the relative values, i.e., the ratio of the heat flux of any model to the 
actual or non- equilibrium flow heat flux, the characteristic under- 
predictive nature of the frozen flow model (ratio = 0.634) and the 
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Table 19. Comparison of the Surface Heat Flux Values for Each Flow Model and Two 

Pyrolysis Gas Compositions at a Front Surface Temperature of 1650°K. 

Pyrolysis Gas Composition 

Mole% 

Flow Model 

Surface Heat Flux 
KJ / M 2 -Sec 

Relative Heat 
Flux, 0 S 

Hydrogen 

Methane 

Carbon Monoxide 

Carbon Dioxide 

Nitrogen 

Water 

Phenol 

31.7 

46.3 

5.5 

2.5 
14.0 

0.0 

0.0 

Frozen 

Non-Equilibrii 

Equilibrium 

574.20 
im 1184.04 
2983.75 

0.485 

1.000 

2.520 

Hydrogen 1 

Methane 

Carbon Monoxide 

Carbon Dioxide 

Nitrogen 

Water 

Phenol 

33.4 

6.7 

3.7 
1.1 
0.0 

48 <9 
6^2 

Frozen 

Non-Equilibrii 

4 

Equilibrium 

771.65 
im 1217.48 

2118.71 

0.634 

1.000 

1.720 

Note: The relative heat flux is defined 

model to the heat flux calculated 
closeness of this value, 0 S , to 1 
(frozen or equilibrium) predicts 
underprediction, while values abo 1 

as the ratio of 
by the non-equi 
determines the 
the actual behav 
ve one indicate 

the heat flux calcu 
Librium flow model, 
iloseness with which 
Lor. Values less th 
in overprediction of 

lated by any 
The relative 
the model 
an one show an 
the true state. 










extremely overpredict ive behavior of the equilibrium flow model (ratio = 
1.720) are shown. It further illustrates that although the frozen flow 
model can accurately describe the true behavior in some cases (low temper- 
atures or high mass flux rates), the equilibrium flow model is totally 
inadequate over the temperature (533- 1925 °K) and mass flux (0.00015- 
0.50 kg/irr-sec) values studied in this research. Therefore, in order 
to obtain an accurate prediction of the energy transfer, a non- 
equilibrium flow model must be used within the transition region. This 
is only possible by considering chemical reactions between the species to 
occur at finite reaction rates as described by reliable kinetic data. 

Again the importance and application of the non- equilibrium flow model 
has been demonstrated and the limitations of the two ideal models shown. 

This discussion forms the basis for extending this research to higher 
temperatures involving more complex reactions and increased number of species. 


Parametric Study of the Flow of Pyrolysis Gases in the Char Zone 

A comparison of the non- equilibrium flow results with the experimental 
data was important in determining the accuracy of the flow model. How- 
ever, very little quantitative information, beyond the discrete sets of 
data for each experiment, was assembled regarding the effect of changing 
mass flux and/or temperature. As a result, a parametric study was under- 
taken to accurately relate the changes in these variables with variations 
in energy absorption within the char. To do this the -non-iterative 
TEMPBE System (NIT) was used in which the back surface temperature and 
temperature gradient were specified as boundary conditions for various 
values of the mass flux. The results of the' calculation were in the 
form of the net heat transfer at the surface, called the approximate 
aerodynamic heating, which was the sum of the surface heat flux and 
radiant heat flux resulting from the calculated front surface temperature: 

q a = -k e dT| +.«otJ (49) 

dz I 

z=L 


where a value of 0.95 was used for the emissivity. This information is 
shown in Figure 10 in which the mass flux is plotted against the aerodynamic 
heating for various heats of pyrolysis, qp, (function of the temperature 
and gradient at the back surface): 


q 


P 



dT I 

dz I z »0 


(50) 
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(Multiple Aerodynamic Heating in BTU/Ft -Sec by 11 to obtain KJ/m -Sec) 



0.02 0.03 

Mass Flux (lb/ft^-sec) 


0.10 0.15 

Mass Flux (Kg/M 2 -Sec) 








where qp is the sum of the energy absorbed by the decomposition of the 
polymer and the energy conducted through the virgin plastic. Results for 
the frozen, equilibrium and non- equilibrium flow models are presented. 

This form of presenting the results is a very convenient and informative 
method as will be seen. 

In a reentry problem one of the important questions asked is what 
is the required heat shield weight for protection for a certain mission. 
Specification of the type of heat shield material to be used (e.g., nylon- 
phenolic resin) brackets the heat of pyrolysis value, while the trajectory 
calculations determine the amount of aerodynamic heating that can be 
expected. For example, an approximate aerodynamic heating rate of 500 
BTU/ft^-sec (5x10^ J/m?-sec) and a heat of pyrolysis of 400 BTU/lb 
(8xl05j/kg) locates three distinctive points on Figure 9; one for each of 
the frozen, equilibrium and non- equilibrium flow models. This corresponds 
to three distinctive values of the mass flux; 0.017 lb/ft2-sec 
(0.085 kg/m?- sec) for frozen, 0.002 lb/ ft?- sec (0.01 kg/m^-sec) for 
equilibrium and 0.009 lb/ft^-sec (0.045. lb/m?~ sec) for non- equilibrium. 

The non- equilibrium flow model accurately predicted the behavior and 
would specify the exact heat shield weight (function of the mass flux) 
required. The frozen flow model shows an over- predict ion because important 
endothermic reactions were omitted, and the equilibrium flow model shows 
an under- predict ion, because reactions were assumed to occur at a greater 
extent than the actual behavior. 

The results presented in Figure 10 also provide a way of determining 
at what point the non- equilibrium flow model changes from the frozen flow 
behavior to a truly non-equilibrium flow condition governed by finite 
reaction rates (6^1(Aj/m2- sec) . This figure very graphically illustrates 
the differences in each model and permits the presentation of a large 
volume of information in a clear and readily accessible manner. 


Calculation of the Reacting Gas Heat Capacity 


In addition to the above information, the reacting gas heat capacity 
for the non- equilibrium flow of pyrolysis products through the char has 
been determined also. This term is very useful in the calculation of 
the one-dimensional, transient response of an ablative composite. The 
energy equation for the transient case can be put in the following form 
(18) for the char zone. 


<2-) c p + 


K+l 

t-i E i s i 

«o(g) 
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Sz Sz 


r ll 

p St 


(51) 



where W is the mass flux of pyrolysis products at z and W D is the mass 
flux of pyrolysis products entering the char. The term in brackets 
is referred to as the effective reacting gas heat capacity. Hence, the 
flow within the char zone can be considered frozen (EH^R^ i= 0) by 
introducing the reacting gas heat capacity as an input function to the 
transient calculations. In Figure 11 a plot of the reacting gas heat 
capacity as a function of temperature is shown for frozen, equilibrium 
and non-equilibrium flow within the char layer up to 3000°F. (1925°K) . 
These curves were calculated for a mass flux of 0.05 Ib/ft^-sec (0.25 
kg/m2-sec), a back surface temperature of 500°F (533°K), and, char 
porosity and thickness of 0.8 and 0.25 inches (0.0064 m) , respectively. 
The differences in the_manner used to calculate the energy transfer by 
chemical reaction (EH^R^) for equilibrium and non- equilibrium flow 
causes the curves to separate as shown. 


Flow of Pyrolysis Products Through Porous Graphite 

There are two important reasons for using porous graphite to 
simulate low density, nylon- phenolic resin chars used in ablative heat 
shield applications. These are availability and machinability of the 
graphite. 

Nylon- phenolic resin chars were obtained from the National Aeronautics 
and Space Administration's Langley Research Center for use in the Char 
Zone Thermal Environment Simulator. Electric air arc jets were used to 
char the nylon- phenolic resins, and this represented a considerable effort 
in supplying just a few specimens for use in this research. Although 
two sections of char were usually obtained from each specimen, the demand 
for additional chars could not be met. In addition to the problem of 
obtaining the samples, the brittle nature of the chars caused serious 
problems in mounting on the char holder section. These complications 
led to the testing and use of graphite, especially, for use in radio- 
active tracer and catalyst evaluation studies. However, for the 
successful substitution of graphite for the chars to be complete, the 
chemical behavior of the two must be essentially the same. This is to 
say that the same reactions and kinetic data important in predicting 
the energy transport within porous chars, must also do the same for 
energy transport in porous graphite. This will be shown in the following 
sections by comparing the exit gas compositions from the Char Zone 
Thermal Environment Simulator with the composition predicted by the 
non- equilibrium flow model, and, by direct comparison of char and 
graphite experimental results over a range of mass flux rates and front 
surface temperature common to both systems. 

Comparison of reacting flow through chars and graphite . -Other than 
the differences in the structural properties of chars and graphite, the 
one most important consideration that must be accounted for is the change 
in mass flux rates caused by differences in material porosity. Therefore, 
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Figure 11- Plot of the Reacting Gas Heat Capacity for Frozen, Equilibrium 
and Non-Equilibrium Flow of Pyrolysis Gases Through the Char Zone. 



to put the materials on a common basis for discussion, the mass flux rates 
previously discussed in terms of the total area (kg/nr total” sec ) must be 
divided by the porosity to obtain rates within the pores (kg/m 2 vo i<ig“Sec) . 
Even though the superficial mass flux rates are different for each porous 
medium, the mass flux rate within the pores will be the same. Again, the 
porosity of the chars and graphite were 0.8 and 0.5, respectively. 

The foregoing discussion is not intended to mean that chemical 
reactions will not be influenced by other structural properties beside 
the porosity (e.g., crystallicity , permeability, etc.). However, it 
emphasizes that the mass flux within the pore spaces must be equivalent 
for a valid comparison. Differences in the chemical reaction rates 
resulting from differences in structural makeup could eliminate graphite 
as a suitable substitute. However, this could only be determined by 
a comparison of experiments conducted over the range of conditions 
for which the proposed flow model is valid. 

In Table 20 the exit gas composition from the Char Zone Thermal 
Environment Simulator for the flow of pyrolysis products through 
graphite are presented for mass flux rates of 0.0034 to 0.0059 lb/ft -sec 
(0.017-0. 030 kg/m 2 -sec) at a front, surface temperature of 
approximately 1950°F (1340°K) . As in the case with chars, there is 
a significant amount of chemical reactions occurring in the porous 
medium for the lower mass flux (0.05 kg/m 2 -sec). More importantly, 
however, is the agreement within experiment error between the measured 
exit gas compositions and the predicted values by the non-equilibrium 
flow model using the same kinetic data employed for the char experiments. 

A similar comparison is presented in Table 22 for an average mass 
flux rate of 0.0035 lb/ft^-sec (0.017 kg/m^-sec) and front surface 
temperatures of 1950 (1340°K) and 2065°F (1410°K) . Again, excellent 
agreement within the experimental accuracy of the analyses was ob- 
tained between the non-equilibrium flow model compositions and the 
experimental values. However, a closer inspection of the graphite 
experiments in Table A-2 of Appendix A shows that the non- equilibrium 
predicted compositions do not agree in some cases with the experimental 
data. This is especially true for experiments XXI-66, XXIII-72, XXIII-73, 
and XXIV-76. These differences were a result of experimental difficulties 
associated with maintaining a steady liquid (water and phenol) feed 
rate as indicated in the experimental summary sheets of Table A-l; and, 
therefore, should not be interpreted as a failure of the model to 
predict the flow behavior in porous graphite. 

In addition to these irregularities in the water-phenol-gas 
experiments, similar disagreement was also observed in several water- 
phenol free pyrolysis product investigations. In experiments IX and XI 
shown in Table 22, increases in methane and corresponding decreases 
in hydrogen concentrations clearly deviated from the predicted frozen 
flow behavior. These observations were attributed to problems associated 
with the fabrication of the graphite specimens. In the earlier experiments, 
the dust- like graphite "fines" which resulted from the milling process 
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Table 22. Flow of Pyrolysis Gases Through Porous Graphite. Effect of 

Having Graphite "Fines" in the Pore Spaces. 


RUN NUMBER 
MASS FLUX FLOW 

FRONT TEMP MODEL 

BACK TEMP 


FROZEN 
EQUILIBRIUM 
NON- EQUILIBRIUM 
EXPERIMENTAL 


FROZEN 
EQUILIBRIUM 

1020 °K I NON- EQUILIBRIUM 
700 °K ! EXPERIMENTAL I 24.6 


IX- 2 7 FROZEN 

0.2850 EQUILIBRIUM 

1000 °K NON- EQUILIBRIUM] 35.0 
650 °K EXPERIMENTAL 27.0 


FROZEN 35.0 

0.0045 EQUILIBRIUM 76.3 

995 °K NON- EQUILIBRIUM 35.0 

705 °K EXPERIMENTAL 28.5 


XI-32 FROZEN 

0.0195 EQUILIBRIUM 

1025 NON-EQUILIBRIUM 

740 °K EXPERIMENTAL 


XI-33 FROZEN 

0,2850 EQUILIBRIUM 

1005 °K NON-EQUILIBRIUM 
. 570 °K EXPERIMENTAL 


43.6 

55.2 


43.6 

4.9 

43.6 

52.2 


w 2 

MOLE MOLE MOLE 


MOLE 

7. 


0.0 


C 6 H 6° C 2 H 4 
MOLE MOLE 


MODEL 


0.0 2.5 
0.0 2.5 
0.0 2.5 
0.0 4.5 


10.75 

84.50 

10.75 


4.51 

52.69 

4.51 



































































were not removed from the pore spaces within the plug. These "fines" 
provided additional surface area which made them highly susceptible to 
reaction with hydrogen to form methane. After complete reaction of the 
"fines" (about five to ten minutes), the hydrogen and methane concentrations 
in the exit stream returned to their initial compositions indicating 
frozen flow behavior. This is illustrated graphically in Figure 12 
in which hyrodgen and methane concentrations (in mole percent) are plotted 
against experimental run time (in minutes) for experiment IX. 

Methane production of this kind was eliminated in subsequent 
experiments (XII, XIII, and XV in Table 23) by first passing nitrogen 
or helium through the graphite plugs prior to mounting in the char 
holder. This removed the "fines" from the pore spaces and eliminated 
the sharp increase in methane observed in Figure 12. The contrasting 
concentration profiles as a function of run time are shown in Figure 13 
for experiment XIII. 

As a final comparison, several char and graphite experiments are 
examined in Table 24. To accomplish this the mass flux rates based on 
the void area were calculated and are shown in brackets. The listing 
in Table 24 is, also, made in order of increasing chemical reactions; 
i.e., low temperatures and high mass flux rates appear first. As seen, 
the char and graphite experiments are indeed compatible and behave, 
from the overall chemical viewpoint, as one material. 

This permits the use of the more easily workable graphite materials in 
experiments designed to study carbon deposition and product distribution 
using Carbon- 14 tracers, and, to investigate the effect of catalysts in 
accelerating the rates of chemical reactions within the char layer. 

Details of the results from these additional experimental investigations 
follow. 


Radioactive Tracer Studies Using Porous Graphite 

Radioactive methane and phenol were used in separate experiments to 
determine the specific products of decomposition from each labeled species. 
Also, the amount and location of carbon deposition in the char due to the 
thermal cracking of each species was determined. 

The method used involved the sampling of the exit gas stream followed 
by gas chromatographic analysis. The fractionated gas chromatographic 
effluent stream was then passed through a combustion furnace forming 
carbon dioxide and water. After trapping the water, the carbon dioxide 
was absorbed in a one molar hydroxide of hyamine (in methanol) solution. 
Collecting the carbon dioxide over small intervals (one half to one minute) 
produced radioactive concentrations corresponding to the separated gases 
indicated on the gas chromatogram. By comparing the two curves for 
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Figure 12. Typical Methane-Hydrogen Concentration Profile in the Exit Gas Stream 
from the Char Zone Thermal Environment Simulator for the Inclusion of 
Graphite "Fines" in the Pore Spaces. 





Table 23 Flow of Pyrolysis Gases Through Porous Graphite. Exit Gas Compositions 
Resulting From Experiments In Which The "Fines" Were Cleared From The 
Pore Spaces By Passing Helium Through the Specimens. 


31 


RUN NUMBE 
MASS FLUX 
FRONT TEMP 
BACK TEMP 


XIII-38 

0.0045 

1195°K 

930°K 


XIII-39 

0.0195 

1210°K 

940°K 


XI II -40 
0.0450 
1230°K 
920°K 


905°K 


FLOW 

MODEL 


FROZEN 
EQUILIBRIUM 
NON- EQUILIBRIUM 
EXPERIMENTAL 1 36.1 


FROZEN 
EQUILIBRIUM 
NON- EQUILIBRIUM 
EXPERIMENTAL 


FROZEN 
EQUILIBRIUM 
NON- EQUILIBRIUM 
EXPERIMENTAL 


FROZEN 31.7 
EQUILIBRIUM 82.3 
NON-EQUILIBRIUM 31.7 
EXPERIMENTAL 30.8 


FROZEN 
EQUILIBRIUM 
NON-EQUILIBRIUM 
EXPERIMENTAL |31.7 


FROZEN 
EQUILIBRIUM 
NON -EQUILIBRIUM 
EXPERIMENTAL 132.5 


'■’“4 

MOLE 

% 


43,6 

42.9 


MM 


M IIEqie&i 

■I 

\ vmrnwnm \ 


46.3 

45.8 


46.3 

46.8 
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MOLE 

% 


MOLE 
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.0 
.6 
.0 
0.0 


.0‘ 
.2 
.0 
0.0 


0.0 

0.1 

0.0 

0.0 


0.0 

0.1 

0.0 

0.0 


0.0 

0.1 

0.0 

0.0 


0. 
0. 
0. 
0.0 


MODEL 


MOLE 
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0.0 

0.0 

0.0 

0.0 


GEH3Q1EQKB 

QiliSp IPtSI 


3.0 

5.2 

4.5 
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Inlet Gas Composition 


Experiment XIII (Graphite) 
Feed Mix 2 
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Figure. 13. .Typical Methane-Hydrogen Concentration Profiles in the Exit Gas Stream 
from the Char Zone Thermal Environment Simulator after Clearing the Pore 
Spaces of Graphite '’Fines’' with Helium or Nitrogen Flow. 




Table 24 


Flow of Pyrolysis Gases Through Graphite and Chars. Comparison of Results 


i 

Ul 

"j 

■ 


RUN NUMBER 
MASS FLUX: 
W 

(w p ) 

FRONT TEMP 
BACK TEMP 

FLOW 

MODEL 

H 2 

MOLE 

7. 

CH. 

4 

MOLE 

7. 

CO 

MOLE 

7. 

C0 2 

MOLE 

7. 

N 2 

MOLE 

% 

H 2° 

MOLE 

% 

W> 

MOLE 

% 

C 2 H 4 

MOLE 

% 

C 2 H 2 

MOLE 

°L 

C-XVITI-56 

FROZEN 

28.9 

6.4 

3.3 

0.8 

0.0 

53.3 

6.9 

0.0 

0.0 

0.01040 

EQUILIBRIUM 

59.7 

3.3 

26.7 

4.0 

0.0 

6.3 

0.0 

0.0 

0.0 

(0.01300) 

NON- EQUILIBRIUM 

23.0 

11.1 

2.8 

2.1 


53.8 


0.0 

0.0 

1190°K 

EXPERIMENTAL 

24.8 

10.3 

3.4 

1.3 • 

0.0 

53.7 


■ tHtH 

- 

820°K 










■i 

G-XXVIII-93 

FROZEN 

25.8 

MEM 

4.3 

1EM 

0.0 

50.1 

11.3 

0.0 


0.0170 

EQUILIBRIUM 

62.7 


33.5 

EM 

0.0 

1.7 


0.0 

0.0 

(0.0340) 

NON -EQUILIBRIUM 

28.8 

10.0 

11.8 

mam 

0.0 

32.6 

10.2 

0.0 

0.3 

1340°K 

EXPERIMENTAL 

32.2 

11.2 

9.2 


0.0 

33.2 

9.8 

0.2 

■ilM 

865 °K 





■■ 






C-XIX-60 

FROZEN 

29.3 

6.0 

3.3 

0.9 

0.0 

53.7 

6.8 


0.0 

0.01240 

EQUILIBRIUM 

62.7 

0.8 

35.2 

1.5 

EHiX 

0.9 

0.0 

0.0 

0.0 

(0.01550) 

NON-EQUILIBRIUM 

38.2 

17.2 

25.2 

8.8 

0.0 

4.9 

5.1 


0.5 

1370°K 

EXPERIMENTAL 

41.3 

18.2 

23.4 

7.6 

0.0 

5.0 

4.5 



1020°K 











G-XXIII-71 

FROZEN 


4.9 

2.9 

0.9 

0.0 

54.7 

10.8 

0.0 

6.0 

0.01550 

EQUILIBRIUM 


0.7 

35.8 

0.3 

0.0 

0.7 

0.0 

0.0 

0.0 

(0.03100) 

NON-EQUILIBRIUM 

41.4 

8.4 

23.9 

5.1 

0.0 

12.5 

8.2 

0.0 

0.6 

1380°K 

EXPERIMENTAL 

39.3 

8.8 

27.2 

5.9 

0.0 

11.2 

7.6 

0.2 

0.3 

865°K 











C-XX-bJ 

FROZEN 

28.6 

5.5 

3.4 

1.0 


53.8 

7.7 

0.0 

0.0 

0.01120 

EQUILIBRIUM 

63.1 

0.3 

36.4 

0.1 

0.0 

0.2 

0.0 

0.0 

0.0 

(0.01400) 

NON-EQUILIBRIUM 

30.7 

21.4 

38.6 

3.4 

0.0 

0.0 

4.6 

0.3 

1.0 

1535°K 
905 “K 

EXPER IMANTAL 

30.4 

19.8 

38.0 

5.7 

0.0 

0.0 

5.0 

0.5 

0.6 


MODEL 


|AP 


cz 


2.0 

2.0 

2.0 

4.0 


10.2 
174.2 
10 . 9 | 


1.0 

1.0 

1.0 

2.0 


21.9 

312.4 

25.6 


3.5 

3.0 

3.0 

3.0 


12.9 

143.0 

26.6 


1.0 

1.0 

1.0 

3.5 


22.0 

350.4 

55.5 


3 . 0| 19.7 
2.51565.0 
2.5! 34.9 
3. Of 































































identical retention times, the relative amount of each carbon- containing 
species formed from the thermal degradation of the labeled pyrolysis 
product entering the char was determined. 

Typical results for Carbon- 14 labeled methane are shown in Figure 
14, in which the gas chromatogram (s) and corresponding radioactivity curve 
are presented. The particular results are for experiment XXIX in which 
the front surface temperature was 1935°F (1333°K) and the gas mass flux 
was 0 . 00591b/ ft 2- sec (0.030 kg/ m 2 -sec). By comparing the two curves, 
the products of methane decomposition were found to be unreacted methane, 
carbon monoxide, carbon dioxide, ethylene and acetylene. These results 
for methane are very important in the light of predicting the manner in 
which energy can be absorbed by chemical reaction. Ethylene and acetylene, 
for example, are indirect products of methane decomposition predicted by 
reactions (39) through (42) in Table 5 , while carbon monoxide and dioxide 
are formed by the reaction of steam with deposited carbon in reactions (46) 
through (48). This information established that the chemical reactions 
used to predict the phenomena occurring in any system are correct. 

A similar discussion is presented for labeled phenol. These results 
are likewise shown in Figure 14. Conditions for the presented data were 
a front surface temperature of 1960°F (1350°K) and a mass flux rate of 
0.0034 lb/ ft 2- sec (0.017 kg/m2-sec). The exit gas products for phenol 
degradation are methane, carbon monoxide, carbon dioxide, ethylene and 
acetylene, as well as unreacted phenol analysed in the liquid phase. 

Once again insight into the kind of reactions necessary to produce the 
products was obtained. The formation of hydrogen and carbon by 
reaction (44) and (45) is probable by the observed carbon deposition 
within the graphite. Hydrogenation of carbon by reaction (43) to form 
methane, followed by the steam-gas reactions (46,47,48) and the hydro- 
carbon cracking reactions (39, 40, 41, 42) accounts for each radioactive 
species observed. 

In both methane and phenol degradation, thermal decomposition of 
the major species in the simulated pyrolysis product stream was described 
and accounted for by the reactions considered important between 500-3000°F 
(533-1925°K) . 

In addition to the product distribution resulting from the thermal 
degradation of methane and phenol, deposited carbon was also observed to 
occur. The location of the carbon deposition within the char layer is 
important in defining the temperature at which reactions become significant. 
This topic is discussed in detail in the next section. 


Carbon Deposition Studies by Radioactive Tracer Methods 

The location and extent of carbon deposition resulting from methane 
and phenol decomposition was determined using Carbon- 14 labeled methane 
and phenol. In the specific cases studied, labeled methane and phenol were 
fed separately as components in the simulated pyrolysis product stream 
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Figure 14. Gas Chromatograms and Radioactivity Curves of the Exit Product Stream from 
the Char Zone Thermal Environment Simulator; (a) C^ labeled methane, (b) phenol. 




entering the char. The char was removed after each experiment and 
sectioned by removing thin layers using emery paper. These layers 
varied between one and ten percent (by weight) of the total char and 
were combusted separately with collection of the carbon dioxide in one 
molar hydroxide hyamine (in methanol) solution. The radioactivity of 
each thin layer was determined and plotted as a function of char depth. 

In Figure 15 such a curve is shown for the thermal degradation of phenol 
and Figure 16 is a similar curve for the decomposition of radioactive 
labeled methane. The hashed- in rectangular blocks represent the total 
percent radioactivity of the thickness of the individual slices analysed, 
while the dotted curve represents the percent radioactivity per unit 
thickness at a particular char depth. The results in Figure 15 are for 
phenol decomposition at a mass flux rate of 0.0059 l'b/ ft 2- sec (0.030 
kg/m^-sec) and a front surface temperature of 1960°F (1350°K). Depo- 
sition of carbon appears to start at a char depth of 0.38 corresponding 
to a temperature of 1300°F (978°K), and continues uniformly to 0.925 
where the temperature is 1925°F (1325°K). At this point a rapid decrease 
is noted indicating either no further carbon deposition or disappearance 
of carbon by chemical reaction. 

Similar results are observed for carbon deposition by methane de- 
composition in Figure 16. Since carbon deposition by methane and/or 
phenol degradation is an increasing function of temperature, and, since 
a substantial amount of phenol and methane is present in the exit gas 
stream, it is unlikely that carbon deposition reactions have terminated. 
Instead, the reaction of the deposited carbon with steam (or carbon 
dioxide) is a more probable explanation of the decline noted in Figure 
16. This is also substantiated by the rapid decrease in water concentration 
at the same temperature where carbon deposition declines. Additionally, 
carbon was observed on the quartz cover plate and inside surfaces of the 
outer char holder section which indicated that the carbon deposition 
reactions were continuing after the gases had left the char surface. 
Therefore, a very comprehensive picture of carbon deposition with 
regard to its location, the causes for its appearance and disappearance, 
and its effect on the exit gas product distribution was obtained. This 
was one additional, important use of an accurate mathematical model in 
predicting phenomena very difficult and often impossible to determine 
with experimental techniques only. 

The combination of the radioactive tracer techniques and the 
non-equilibrium flow analysis will be applied in evaluating various 
catalysts for accelerating the chemical reactions and, thereby, increasing 
the energy absorbed within the char zone. The effectiveness of each catalyst 
will be determined in the following section by comparing the results 
with data from non-catalytic experiments. 
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Catalytic Reactions of the Pyrolysis Products in the Char Zone 

The discussion thus far has shown that chemical reactions within the 
char layer are very important modes of energy absorption. It was also 
pointed out that chemical non-equilibiium becomes important between 
2000 - 2500 °F (1367 - 1644 e K) for a mass flux rate of 0.05 lb/ft 2 -sec 
(0.25 kg/m 2 - sec) (Tables 15 through 18). Below this range, the flow of 
pyrolysis products through the char is essentially frozen, while above 
2000 n F (1367°K), the flow is best described by finite reaction rates for 
the important chemical reactions taking place within the char between the 
pyrolysis products. 

In reentry applications where the temperature gradient across the 
ablator may vary from 500 - 5000°F (533 - 3033°K), the frozen state (less 
than 1367°K) can extend over nearly one third the entire thickness. 

Within this region heat absorption by conduction and convection are 
the only important modes of energy transfer. In order to increase the 
capacity of this region to absorb heat, and thus proportionately reduce 
the total heat shield weight requirement, the introduction of a catalyst 
to initiate chemical reactions in the lower temperature range (<1367°K) 
was evaluated. 

There are two types of catalyst systems: homogeneous and heterogeneous. 

The first involves the introduction of a chemical compound which is in 
the same phase as the pyrolysis product stream. The homogeneous catalyst 
effectively reduces the energy of activation by interacting with the 
various species present. To illustrate this interaction the catalytic 
chlorination of nitrous oxide in the presence of bromine is used as an 
example ( 74 ) of the action of a homogeneous catalyst. The direct 
chlorination occurs by reaction (52): 

2NO + Cl 2 -♦ 2NOC1 (52) 


Introducing bromine results in a two reaction sequence as follows: 


2N0 + Br 2 -♦ 2N0Br 

(53) 

2N0Br + Cl 2 -+ 2N0C1 + Br 2 

(54) 


The reaction with the greatest activation energy between (53) and (54) 
is the rate determining step for the sequence. Bromine is considered a 
catalyst if, and only if, the energy of activation of the rate determining 
step is smaller than the energy of activation of reaction (52). Such is 
the case for this example. Other homogeneous catalysts are iodine, hydrogen 
bromide, hydrogen fluoride, nitric oxide, chlorine and mercury ( 75 , 76, 77 , 
78, 79). 
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In the case of the halogens, incorporation of an organic halide into 
the composite which thermally degrades at or near the temperature for the 
nylon-phenolic resin could be used to introduce the catalyst into the 
pyrolysis product stream. One example for bromine is tribromobutane which 
vaporizes at 500°K (nylon-phenolic resin at 525°K) and forms HBr and 
cracked products of an olefin. 

The second kind of catalyst system with application to flow of 
hydrocarbon products through porous media are the heterogeneous 
catalysts. These consist of a thin dispersion of an active metal on 
a porous solid, called a catalyst support. Gases are thus absorbed on 
the metal surfaces as they flow through the porous solid, undergo chemical 
reactions, and desorbe back into the gas stream. The kind of solid 
supports used vary from clays and alumina to porous carbon. 

Heterogeneous catalysts are widely used in the petroleum and 
chemical industries for accelerating hydrogenation and dehydrogenation 
reactions, hydrocracking reactions and hydroforming reactions. Some 
typical active metals used in these applications are platinum, tungsten, 
molybdenum, palladium, etc. (80). One example involves the catalytic 
hydrogenation of benzene to cyclohexane at room temperature with platinum 
on porous carbon supports (80). Cyclohexane is then cracked to lower 
molecular weight compounds at 722°K (81). Details for the catalytic 
cracking of numerous organic compounds are presented by Vogh (82). 

In many cases the use of heterogenous catalysts is restricted 
to applications which do not contain compounds that deactivate the 
metal surfaces. Some of these so called poisons are carbon monoxide, 
sulfur and deposits of carbon or coke. Although there is no sulfur in 
the pyrolysis product stream, carbon monoxide and deposited carbon are 
present requiring additional screening of the heterogeneous catalyst 
considered. The activity and selectivity of the catalyst chosen for this 
application is discussed in a subsequent section. 

Because the heterogeneous catalyst exists in a different state, 
the solid phase, then the pyrolysis products, introduction into the system 
requires more detailed planning. There are two possible techniques that 
permit the placement of a finely dispersed metal catalyst on the char 
of a charring ablator. The first takes advantage of a nylon- platinum 
catalyst used to hydrogenate benzene to cyclohexane (83). The inclusion 
of this platinum impregnated nylon with nylon and phenolic resin during 
the molding process could be made. During ablation, the nylon would 
degrade and release the metal catalyst which would be distributed on the 
surface of the formed char layer. The presence of water and hydrogen at 
the lower temperatures (<750°K) would prevent coking, leaving the metal 
sites exposed to promote the pyrolysis reactions. 
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The second method is similar to the method used for introducing a 
homogeneous catalyst into the pyrolysis product stream- In this case an 
organo-metallic compound such as nickel stearate (84) could be included 
which would vaporize in the decomposition zone with deposition of nickel 
on. the char surface. This action is commonly observed in vapor phase 
cracking processes (84,85) in which increased activity of the cracking 
catalyst results in excessive carbon and hydrogen formation. Other 
similar compounds are the carbonyl compounds of nickel, iron and cobalt 
(86) . The combination of both catalysts systems may also be possible 
by using compounds containing both metal and halogen atoms, such as 
platinum iodide. The advantage of this type of co-catalyst would exist 
only if both groups were found to accelerate chemical reactions within 
the char layer. 

The following sections will describe the results obtained in tests 
using the Char Zone Thermal Environment Simulator. Each catalyst system 
will be compared with non-catalytic data to determine the extent of 
chemical reaction due to the addition of the catalyst. The preparation 
of each catalyst and the procedure for introducing it into the experimental 
simulation will also be discussed. 

Homogeneous catalysis of the pyrolysis products . -Uni ike the 
heterogeneous catalysis systems, very little information regarding the 
activity of various homogeneous catalyst systems is contained in the 
literature. One source, however, reported the relative activity of 
several organic halides and halogen catalysts for the catalytic degradation 
of hydrocarbons to carbon monoxide, carbon dioxide and organic acids, 
aldehydes and ketones (87) . A list of the relative activities of these 
catalysts are presented in Table 25 with iodine as a reference (relative 
activity of 100). No relative activities of nitrous oxide or mercury 
were found. Also hydrogen fluoride was omitted from the above list of 
relative activities. 

Although the above may indeed be excellent homogeneous catalysts, 
certain aspects of the ablative process prohibit their use. For example, 
nitrous oxide, while exhibiting excellent catalytic activity for the 
thermal degradation of certain hydrocarbons requires a concentration 
too great to be practically included in the composite (75) . Similarly, 
hydrogen fluoride and mercury are almost exclusively used as liquid phase 
catalysts (88, 89). As a result, attention in this research was given 
to the halides which were not only reported as good catalysts in hydro- 
carbon decomposition and oxidation reactions (90), but also required in 
small enough concentrations to be conviently and economically feasible 
for ablative heat shield applications. 

Experimental results for homogeneous catalysis of the pyrolysis 
product reactions with bromine . -As a matter of convenience in the 
experimental simulation apparatus, bromine was selected as a representative 
halide catalyst. It was convenient to dissolve bromine in the water and 
feed the resulting solution to the char holder with phenol and the gaseous 
pyrolysis products. The concentration of bromine in the water solution 
was varied from 1.0 to 4.0% by weight. 


- 65 - 



I 

O'* 

o^ 

I 


Table 25. 


Relative Activity of the Halogens and Organic Halides as Catalysts in 
the Thermal Degradation of Pyrolysis Gas Products. 


Conditions : Oxidation of Aldehydhes and Ketones at 400-500° C and one Atmosphere. 


Component Name 


Relative Efficiency 


Iodine (Reference) 
Isopropyl Iodide 
Ethyl Iodide 
Methyl Iodide 
Hydrogen Bromide 
Isopropyl Bromide 
Ethyl Bromide 
Dibromoacetylene 
Chlorides 





The selection of bromine as a suitable and typical catalyst was 
based on evidence of strong aromatic ring attack, such as with phenol, 
to form tribromophenol at ambient conditions. In addition it has ex- 
cellent oxidative properties in decomposing hydrocarbons such as 
propane to carbon monoxide, carbon dioxide and other organic species 
such as organic acids, aldehydes and ketones (90). Furthermore, the 
formation of HBr from hydrogen and bromine occurs at moderate temper- 
atures (475-575°K) by the following mechanism (91). This gives the 
additional catalytic benefit of HBr being present. 

Br 2 -> Br + Br 

Br + H 2 -» HBr + H 

H + Br 2 -* HBr + Br 

H + HBr -* H 2 + Br (55) 

Br + Br -» Br 2 * 

Br 2 * -» Br + Br 

Also it established a reference to the remaining halide catalysts 
contained in Table 25 and thus makes an exhaustive investigation 
unnecessary if significant promotion of chemical reaction rates are 
found. The above mechanism is favored by low pressure and large 
surface area, both of which exist in the char zone during reentry. 

Results of several experiments are shown in Table 26. The 
experimental exit gas compositions are no longer predicted by the 
non- equilibrium! flow model within the experimental accuracy as 
observed for the non- catalytic experiments. Also the reactions are 
not at equilibrium as seen from the table. 

The effect of bromine as a catalyst is better illustrated by 
comparing experiment XXVIII- 92, in which the mass flux was 0.0044 
Ib/ft2-sec (0.022 kg/m^-sec), the front surface temperature of 1920°F 
(1322°K) and no bromine, with experiment XXV- 81 (mass flux of 0.010 
kg/m^-sec, a front surface temperature of 1995°F (1370°K) and 4% (wt) 
bromine catalyst) in Table 27- The conditions are almost the same, 
and the extent of reaction is greater for the experiment with bromine 
present. The addition of bromine accelerated the chemical reactions 
as seen by the facts that significantly more of the water (17 . 6% rather 
than 25.47.) and phenol (7.27 0 rather than 9.57.) had been consumed than 
would have been if bromine had not been present. This is also seen 
by direct comparison of the experimental exit gas compositions from 
each of the similar experiments. A measured decrease in the water 
and phenol concentrations, with corresponding increases in methane, 
carbon monoxide, carbon dioxide, ethylene and acetylene are obtained in 
the bromine catalysed experiment. These same trends are observed in all 
investigations with bromine and are independent of the catalyst 


- 67 - 



Table 26. Flow of Pyrolysis Gases Through Graph: 

the Exit Gas Compositions 


RUN NUMBER 
MASS FLUX 
FRONT TEMP 
BACK TEMP 

FLOW 

MODEL 

EH 

CH. 

4 

MOLE 

% 

CO 

MOLE 

% 

co 2 

MOLE 

% 

XXV-82 

FROZEN 

36.9 

MEM 

4.2 

1.3 

0.0325 

EQUILIBRIUM 

66.2 

■ 

32.1 

0.2 

1380°K 

NON-EQUILIBRIUM 

40.6 

III 

14.1 

2.7 

_86Q°K 

EXPERIMENTAL 

41.6 

HMl 

18.9 

3,4 

XXV 1 1- 87 

FROZEN 

2.3 

3.4 


0.1 

0.0160 

EQUILIBRIUM 

73.7 

1.5 

24.1 

0.3 

1350°K 

NON-EQUILIBRIUM 

2.8 

3.6 

0.5 

0.1 

845 °K 

EXPERIMENTAL 

4.8 

2.7 

0.8 

0.1 

XXVII-89 

FROZEN 

2.1 

3.2 

0.4 

■ 

0.0254 

EQUILIBRIUM 

51.1 

0.6 

46.7 


1365 °K 

NON-EQUILIBRIUM 

18.3 

2.8 

13.3 


850°K 

EXPERIMENTAL 

28.2 

3.4 

21.5 

6.0 

XXVI-84 

FROZEN 

34.3 

6.6 

3.9 

1.1 

0.0065 

EQUILIBRIUM 

65.0 

1.0 

32.5 

0.4 

1370°K 

NON-EQUILIBRIUM 

39.5 

10.7 

16.6 

3.8 

755 °K 

EXPERIMENTAL 

37.8 

15.3 

21.0 

4.8 
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Table 

27. Flow of Pyrolysis Gases Through Graphite. Comparison of Bromine Catalyzed 

and Non-Catalytic Exit Gas Compositions. 

RUN NUMBER 


«2 

m 

CO 






HB 


MASS FLUX 

FLOW 

co 2 

N 2 

H 2 0 

C-H.O 
6 6 

C 2 H 4 


COMMENTS 

FRONT TEMP 
BACK TEMP 

MODEL 

MOLE 

7. 


MOLE 

7. 

MOLE 

7. 

MOLE 

7. 

MOLE 

7. 

MOLE 

7. 

MOLE 

7. 



XXV-81 

FROZEN 

31.5 

6.1 

3.5 

1.1 

0.0 

47.1 

10.8 


0.0 

Bromine 

0.0190 

EQUILIBRIUM 

64.7 

0.9 

33.2 

0.4 

0.0 

0.9 

0.0 

0.0 

0.0 

in Water 

1350°K 

NON-EQUILIBRIUM 

36.4 

9.9 

14.6 

3.7 

0.0 

25.4 

9.5 

0.0 

0.4 

47. 

903°K 

EXPERIMENTAL 

40.2 

11.9 

15.1 

7.0 

0.0 

17.6 

7.2 

0.3 

KH 

(by weight) 

XXVIII-92 

FROZEN 

27.5 

5.3 

3.1 

1.0 

in 

52.5 

10.6 

0.0 

0.0 

No 

0.0220 

EQUILIBRIUM 

62.5 

1.2 

33.9 

0.8 

ESI 

1.6 

0.0 

0.0 

0.0 

Bromine 

1320°K 

NON- EQUILIBRIUM 

29.3 

MEM 

8.6 

2.5 

KS9 

42.2 

10.1 

0.0 

0.2 


855°K 

EXPERIMENTAL 

30.0 

mnm 

8.0 

3.4 

0.0 

39.9 

9.6 

0.2 

0.2 



































concentrations used (1.0 to 4.0% by weight). 

The effect of the bromine catalysis Is also seen by Inspecting the 
carbon deposition profiles for the thermal decomposition of phenol in 
the simulated pyrolysis product stream. This Is presented in Figures 
17 and 18. The carbon deposition profile for the non-catalytic experiment 
is shown in Figure 17 for a mass flux of 0.0059 lb/ft2-sec (0.030 kg/m2-sec) 
and 1935°F (1330°K) , Deposition begins at a distance of 0.38 in the char 
where the temperature is 1312°F (980°K) with a uniform Increase to a 
maximum value at 0.925 or 1915°F (1320°K). At this point the profile 
decreases. In contrast to this curve, Figure 18 represents the carbon 
deposition profile for the bromine catalysed experiment in which the 
mass flux was 0.0075 lb/ ft 2- sec (0.038 kg/m^-sec) and 1900°F (1315°K). 

There is a definite shift in the carbon deposition curve with the 
maximum value moving from 0.925 to 0.71 or from a temperature of 1915 ®F 
(1320°K) to 1,630°F (1160°K). Deposition again starts at a distance of 
0.38 (975°K). Also additional carbon deposition was noted near the 
front surface at a temperature of about 1925°F (1325°K). 

Similar profiles are shown in Figure 19 and 20 for carbon- 14 
labeled methane decomposition. Although the results are not as 
pronounced as the phenol data, a detactable shift in the carbon 
deposition pattern is again observed. Deposition begins at 0.48 (1120°K) 
for the non-catalytic experiment compared with 0.42 (1090°K) for the 
bromine catalyzed case. The point of maximum deposition is shifted 
away from the front surface to a position corresponding to 1893°F 
(1.300°K) or 0.85. The peak for the bromine- free experiment is located 
at 0.98 (1330°K) . 

In summary these results show that bromine is an active catalyst 
for the catalytic cracking of phenol in the presence of hydrogen 
(hydrocracking) and to a lesser extent for the hydrocracking of 
methane. A similar discussion regarding the cracking of the pyrolysis 
products using heterogeneous catalyst is presented in the following 
section. 

Heterogeneous catalysis of the pyrolysis product reactions using 
a tungsten- molybdenum co-catalyst . -A great deal of research into heteroge- 
neous catalysts, their applications and activities has been reported in 
the literature (80,81,82,83). These catalyst are excellent hydrogenation 
accelerators for a number of hydrocarbons common to the petroleum and 
chemical industries. However, they are also susceptible to deactivation 
by reaction, adsorption or coating by several poisons. The two poisons 
which are present in the pyrolysis gases are carbon monoxide and coke 
(or carbon) formation. These two poisons are present in the char zone 
and must be considered when selecting possible heterogeneous catalysts. 

These poisons rule out the use of platinum, paladium, rhodium, nickel and 
selenium since these are all poisoned by carbon monoxide. In light 
of this discussion, tungsten, because of its relatively good activity 
in systems containing carbon monoxide, and molybdenum, because of its 
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Figure 17. Carbon Deposition Profile for the Thermal Decomposition of Phenol, a Major 

Component in the Pyrolysis Gas Stream. 
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Figure 19. Carbon Deposition Profile for the Thermal Degradation of Methane* a Major 

Component in the Pyrolysfs Gas Stream. 
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Figure 20. 




high s©l©ctivlty In the thermal degradation of hydrocarbons, were 
selected as co-catalysts. 

The method used to disperse these metals on the graphite specimens 
was the standard procedure to prepare heterogeneous catalysts and Is as 
follows. First, the metals were obtained as metal acids (anhydrous) 
and added to hot (353°K) sulfuric acid. Molybdenum was completely 
dissolved while tungsten formed a saturated solution. The graphite 
specimens were placed in the hot solution and stirred vigorously for 
thirty minutes. The second phase of the procedure involved the passing 
of carbon disulfide vapors through the graphite to convert the metal 
oxides to sulfides which Increases the catalytic activity of the 
metals. The specimens were then dried at 316 °K and reweighed to 
determine the weight of catalyst dispersed within the pores. The 
catalyst concentration varied from 5 to 6 percent (by weight) of the 
co-catalyst (50:50). This is typical of the dispersed (metal) phase 
composition of heterogeneous catalysts. 

The effect of this catalyst on the reactions of the pyrolysis 
products in the char zone is shown In Table 28. Although there are 
noticeable differences between the experimentally measured exit gas 
compositions and the computed values for the uncatalyzed case, they 
are not as pronounced as the results for bromine. This is better indicated 
by comparing the results of non-catalytic experiments (XXVIII- 92 and XXIX- 
94) with the values obtained in the heterogeneous co-catalyst systems 
(XXXII- 99 and XXXI-98). The exit gas compositions are shown in Table 29 
for an average mass flux rate of 0.0048 lb/ft^-sec (0.025 kg/m^-sec) 
and a front surface temperature range of I860 0 to 1945°F (1290°K- 
1336°K). No detectible difference between the four experiments is 
determined. Similarly, a comparison of the carbon deposition profiles 
for carbon- 14 labeled phenol (XXVIII and XXXII) and methane (XX.IX and 
XXXI) show no substantial shift in the position of the maximum 
deposition or the location where carbon deposition begins. These 
are shown in Figures 21 through 24. A slight difference ir: behavior 
near the front surface is observed. However, this represents no 
substantial change from the non-catalytic behavior. 

A series of experiments with a platinum catalyst dispersed on 
the graphite specimens was also conducted using phenol-water free 
pyrolysis gases having a relatively low carbon monoxide concentration. 

These studies showed no catalytic activity over a temperature range of 
1370° to 1755°F (1015°K to 1230°K). The results of these tests are 
summarized in Table 30. The flow remained frozen in each experiment 
studied . 


Summary of the Catalytic Studies to Increase Reactions in the Char Zone 

Results for adding bromine to the pyrolysis product stream as a 
homogeneous catalyst indicated increased chemical reactions within the char 
zone. This acceleration was shown by comparison of the exit gas composition 
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Table 28 Flow of Pyrolysis Gases Through Graphite. Effect of Molybdenum-Tungsten 

Co-catalyst on the Exit Gas Composition. 


RUN NUMBER 
MASS FLUX 
FRONT TEMP 
BACK TEMP 

FLOW 

MODEL 

«2 

MOLE 

7. 

CH, 

4 

MOLE 

7. 

CO 

MOLE 

7. 

co 2 

MOLE 

7. 

N 2 

MOLE 

7. 

H 2 0 

MOLE 

7. 

C,H,0 
6 6 

MOLE 

7. 

C 2 H 4 

MOLE 

"% 

C 2 H 2 

MOLE 

7o 

XXXI-97 

FROZEN 

20.0 

4.2 

2.1 

0.6 

0.0 

m 

11.0 

0.0 


0.0320 

EQUILIBRIUM 

59.4 

1.7 

33.7 

2.0 

0.0 

mi 

0 t 0 

0.0 

0.0 

1300°K 

NON-EQUILIBRIUM 

20.4 

4.6 

3.0 

0.9 

0.0 

60.1 

10.9 

0.0 

0.0 

740°K 

EXPERIMENTAL 

21.6 

5.6 

4.1 

1.3 

0.0 

58.0 

9.0 

0.2 

0.3 

XXXI-98 

FROZEN 

26.2 


2.9 

0.8 

0.0 

52.7 

IWfli 

0.0 

KSI 

0.0240 

EQUILIBRIUM 

62.6 


33.1 

1.0 

0.0 

1.9 

IQS 

0.0 


1340°K 

NON-EQUILIBRIUM 

28.4 


6.9 

1.5 

0.0 

45.2 


0.0 


770°K 

EXPERIMENTAL - 

29.2 

11.1 

8.7 

3.6 

0.0 

38.0 

ESI 

0.3 

KB 

XXXII-99 

FROZEN 

25.4 

5.4 

2.3 

0.8 

0.0 

53.9 

11.7 

0.0 



EQUILIBRIUM 

62.6 

1.4 

33.1 

1.0 

0.0 

1.9 

0.0 

0.0 

0.0 

1280 K 

NON-EQUILIBRIUM 

26.3 

6.3 

5.0 

1.4 

0.0 

49.6 

11.5 

0.0 

0.1 

730°K 

.. 

EXPERIMENTAL 

27.6 

7.2 

6.6 

2.6 

0.0 

45.2 

10.2 

0.2 

0.4 

XXXII-100. 

FROZEN 

14.0 

2.3 

1.5 

0.4 

0.0 

69.4 

12.4 

0.0 

0.0 

0.0490 

EQUILIBRIUM 

62.6 

1.4 

33.1 

1.0 

0.0 

1.9 


0.0 

0.0 

1330°K 

NON-EQUILIBRIUM 

15.2 

2.4 

2.8 

0.6 

0.0 

66.8 

12.2 

0.0 

0.0 

750°K 

EXPERIMENTAL 

16.6 

3.2 

3.7 

0.8 

0.0 

64.0 

10.8 

0.4 

0.5 


CATALYST 


Molybdenum 
& Tunj^ten 

by weight 
as sulfidesl 


as 

above 


as 

above 


as 

above 
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Table 29 . Flow of Pyrolysis Products Through Porous Graphite. Comparison of the Exit Gas 

Composition for Molybdenum- Tungsten Catalysed and Non-Catalyst Experiments. 


RUN NUMBER 
MASS FLUX 
FRONT TEMP 
BACK TEMP 

FLOW 

MODEL 

m 

m 

CO 

MOLE 

7. 

co 2 

MOLE 

7. 

n 2 

MOLE 

7. 

h 2 o 

MOLE 

7. 

C 6 H 60 

MOLE 

7, 

c 2 h 4 

MOLE 

7. 

c 2 h 2 

MOLE 

7o 

——I 

CATALYST 

XXVIII- 92 

FROZEN 

27.5 

5.3 

3.1 

1.0 

0.0 

52.5 

10.6 

KBS 

0.0 


0.0220 

EQUILIBRIUM 

62.5 

1.2 

33.9 

0.8 

0.0 

1.6 

0.0 

HI 

0.0 

no 

1325°K 

NON- EQUILIBRIUM 

29.3 

MEM 

8.0 

2.5 

0.0 

42.2 

10.1 

HI 

0.2 

catalyst 

855 °K 

EXPERIMENTAL 

30.0 


8.0 

3.4 

0.0 

37.9 

11.6 

KB 

0.2 


XXIX- 94 

FROZEN 

25.6 

4.1 

2.9 

0.8 

0.0 

55.2 

11.5 

0.0 

0.0 


0.0250 

EQUILIBRIUM 

61.8 

1.1 

34.9 

0.7 

0.0 

1.5 

0.0 

0.0 

0.0 

no 

1325 °K 

NON- EQUILIBRIUM 

28.1 

5.8 

7.8 

2.1 

0.0 

45.1 

10.9 

0.0 

0.3 

catalyst 

760°K 

EXPERIMENTAL 

25.5 

5.8 

10.5 

2.7 

0.0 

42.7 

12.2 

0.2 

0.2 


XXXI- 98 

FROZEN 

26.2 

5.5 

2.9 

0.8 

0.0 

52.7 

11.9 

m 

0.0 

RfflfflRHi 

0.0240 

EQUILIBRIUM 

62.6 

1.4 

33.1 

1.0 

0.0 

1.9 

0.0 

139 

0.0 

Tungsten 

1340 °K 

NON- EQUILIBRIUM 

28.4 

6.5 

6.9 

1.5 

0.0 

45.2 

11.4 


0.2 

5-67. (wt) 

770°K 

EXPERIMENTAL 

29.2 

11.1 

8.7 

3.6 

0.0 

38.0 

8.8 

o 

0.3 

as sulfide 


FROZEN 

25.5 

5.3 

2.8 

0.8 

0.0 

53.9 


0.0 

0.0 

Molybdenun 

0,0245 

EQUILIBRIUM 

62.6 

1.4 

33.1 

1.0 

0.0 

1.9 

iBeS«§ 

0.0 

0.0 

Tungsten 

1290°K 

NON- EQUILIBRIUM 

26.3 

6.3 

5.0 

1.4 

0.0 

49.6 

■HH 

0.0 

0.1 

5-67. (wt) 

720°K 

EXPERIMENTAL 

27.6 

7.2 

6.6 

2.6 

0.0 

45.2 

■ 

10.2 

0.2 

0.4 

as sulfide 
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Figure 21 , Carbon Deposition Profile for the Thermal Degradation of Methane, a Major 

Component in the Pyrolysis Gas Stream. 
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Figure 22 . Carbon Deposition Profile for the Thermal Degradation of Methane Catalyzed 

with a Molybdenum-Tungsten Co-Catalyst Coating on the Graphite. 
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Figure 24. 


Carbon Deposition Profile for the Thermal Degradation of Phenol Catalyzed 
with a Molybdenum-Tungsten Co-Catalyst Coating on the Graphite. 



Table 30. Flow of Pyrolysis Products Through Porous Graphite. Exit Gas Compositions 

of Platinum Catalysed Experiments. 


RUN NUMBER 
MASS FLUX 

FLOW 

H 2 

CH. 

4 

CO 

C0 2 

N 2 

h 2 o 

C 6 H 6° 

C 2«4 

c 2 h 2 

FRONT TEMP 

MODEL 

MOLE 

MOLE 

MOLE 

MOLE 

MOLE 

MOLE 

MOLE 

MOLE 

MOLE 


BACK TEMP 


% 

% 

% 

% 

% 

7. 

% 

7 . 

% 

XIII-38 

FROZEN 

31.7 

EBI 

5.5 

2.5 


0.0 

0.0 

0.0 

0.0 

0.0045 * 

EQUILIBRIUM 

82.3 

■oil 

7.0 

0.0 

9.5 

0.2 

0.0 

0.0 

0.0 

1175°K 

NON- EQUILIBRIUM 

31.7 

46.3 

5.5 

2.5 

14.0 

0.0 

0.0 

0.0 

0.0 

924° K 

EXPERIMENTAL 

32.6 

41.3 

8.9 

5.2 

14.0 

0.0 

0.0 

0.0 

m 

XIII- 39 

FROZEN 

31.7 

46.3 

5.5 

2.5 

M 

KH 

B1 

0.0 

0.0 

0.0195 

EQUILIBRIUM 

82.4 

1.0 

7.0 

0.0 

RSI 

0.1 

■ 

0.0 

0.0 

1215°K 

NON- EQUILIBRIUM 

31.7 

46.3 

5.5 

2.5 


0.0 


0.0 

0.0 

930°K 

EXPERIMENTAL 

30.8 

45.8 

6.0 

4.4 

14.0 

0.0 

E89 

0.0 

0.0 

VIII- 22 

FROZEN 

35.0 

43.7 

4.9 

2.9 


0.0 

0.0 

in 

0.0 

0.0190 

EQUILIBRIUM 

79.8 

3.0 

6.7 

0.1 

KU 

0.7 

0.0 

in 

0.0 

1065°K 

NON- EQUILIBRIUM 

35.0 

43./ 

4.9 

2.9 

13.6 

0.0 

0.0 

B9 

0.0 

660°K 

EXPERIMENTAL 

36.2 

41.4 

5.2 

4.2 

13.0 

0.0 

0.0 

0.0 

0.0 










































based on the non-equilibrium (non-catalysed) flow model. Furthermore, 
carbon deposition profiles were used to locate the position and temperature 
where deposition resulting from methane and phenol thermal degradation 
occurred. 

The use of bromine (or HBr) as a representative homogeneous 
catalyst made a detailed investigation unnecessary because of the 
relative activity of various halides to hydrogen bromide in Table 
26. By comparison, iodine and organic iodides should have a greater 
influence on the pyrolysis product reactions, while chlorine or the 
organic chlorides should have a lesser influence. 

No detectable change in the exit product composition, and, there- 
fore, in the rates of chemical reaction of the pyrolysis products, 
were observed for tungsten and molybdenum co-catalysts and platinum 
over a temperature range of 1845° to 2300°F (1280°K to 1535°K). 


Oxidative Degradation of Low Density Nylon- Phenolic Resin Chars 
in an Air Stream at Elevated Temperatures 

The oxidative degradation of low density nylon- phenolic resin 
chars is important in reentry applications as a result of the hot, 
shock heated air streams flowing adjacent to the char front surface. 

It is believed that a portion of this air stream penetrates the char 
front surface forming large cavities and irregular shaped, erroded 
areas on the surface of the heat shield during reentry or in models 
tested in arc jets. In order to determine the extent of air oxidation, 
experiments were made in the Char Zone Thermal Environment Simulator 
in which air, at rates corresponding to reentry values, was injected 
into the char at the front surface. These experiments simulated a 
maximum air mass flux, because there is pyrolysis gas flow counter 
to the air during ablation. 

To determine whether air reacted at the surface or in depth, 
samples of the exit gas stream were taken at short intervals to 
determine the oxygen conversion within the char. The result of 
these analyses for two separate experiments are shown in Figure 25. 

At an air mass flux of 0.0057 lb/ft^-sec (0.029 kg/m^-sec) and a 
front surface temperature of 1485°F (1080°K), the maximum conversion 
of oxygen at the back surface was 53% corresponding to an oxygen 
concentration of 10 mole percent leaving the back surface. A back 
surface temperature of (670°K) was measured for the 1/8 inch (0.0032 m) 
thick char. Increasing the air flow rate to 0.035 lb/ft^-sec (0.18 kg/m^- 
sec) and a front surface temperature of 2047°F (1393°K), increased 
the maximum conversion at the rear surface to 81% or a concentration 
of four percent leaving the char back surface. The back surface, 
temperature was measured at (425 °K) . The char thickness was 1/4 inch 
(0.0064 m) . The relatively low back surface temperatures indicate 
that the highly exothermic reaction probably occurs near the front 
surface. Although oxygen conversion at the back surface was never 
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Experiment XVII- 24 

T Front = 1380°K L = 0.00315 M 
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Tgack = 665 °K 
L - 0.0063 M 
W «=> 0.0285 Kg/ m^- sec 


L5 20 25 30 

Experimental Run Time, Minutes 


Figure 25. Conversion of Oxygen as a Function of Experimental Run Time for the Air 

Oxidation of Porous Graphite at High Temperatures. 




complete indicating the presence of oxygen within the pores at all char 
depths the extent of oxidation in depth appears small in comparison to 
the surface oxidation. The three large erroded areas indicate strong sur- 
face attack. The rounded edges near the errosions further support the 
concept that oxidation occurs where there is the least resistance to 
flow. The time for these holes to be formed can be estimated by comparing 
the length of time necessary to obtain near zero oxygen conversion at the 
back surface in Figure 25. 

In addition to the experimental data 5 a modification to the 
isothermal analysis was made by considering temperatures within 
the char linear with char distance. This approximation permits 
the calculation of the oxygen conversion for a given char depth 
(residence time) and front surface temperature. A comparison 
of the calculated and measured maximum conversions is shown in 
Table 31. The conversion of oxygen with char depth is also presented 
in Figure 26 for the linear approximation of the temperature and 
also the isothermal evaluations at the front and back surface 
temperatures. The data used to describe the air oxidation of 
graphite is presented in Table 31. Agreement of the calculated 
and experiment conversions at the back surface is within the 
experimental accuracy of the analyses. 
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Fractional conversion of Methane 



Dimensionless Char Distance 
(z/L) j 1 “ Char Back Surface 


Figure 26. Convers Ion of Oxygen for the Flow of Air Through a 
Low Density Nylon- Phenolic Resin Char at a Mass Flux of 
0.1750 Kg/mr-sec and a Temperature 'of 1060 °K. 
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Table 31. 


Air Oxidation of Low Density Nylon- Phenolic Resin Chars. A Detailed 
Comparison of the Calculated and Measured Oxygen Conversion at 
the Char Back Surface (Exit Surface) . 


Experiment Number: 

X-29 

XVI I- 54 

Conditions: 



Mass Flux, Kg/m2-sec 

0.0285 

0.1750 

Front Temperature, °K 

1060 

1380 

Back Temperature, °K 

665 

425 

Char Thickness , M 

0.0063 

0.00315 

Experiment Run Time, Min. 

35 

25 

Measured Maximum Conversion: (%) 

43 

83 

Calculated Maximum Conversion: 



(1) at the Front Temperature 

96 

100 

(2) at the Back Temperature 

8 

22 

(3) at a Temperature which varied 

45 

86 

linearly with char depth 



(4) at the Average Temperature 

41 

81 




CONCLUSIONS 


Based on the experimental and theoretical results of this 
research the following conclusions are drawn: 


Reacting Flow of Pyrolysis Products in the Char Zone 

1. The reacting flow of pyrolysis products from nylon- phenolic 
resin composites in the char zone is accurately described by a non- 
equilibrium model employing reaction kinetic data of the important 
reactions occurring in the system. For the simulated pyrolysis 
product compositions in Table 6, mass flux values between 0.00003 - 

0.10 l.b/ft2~sec (0.00015 - 0.5 Kg/m2-sec) and a front surface 
temperature range of 1350 - 2300°F (1000°K - 1535°K), the reactions 
and associated kinetic data in Table 5 accurately describe the 
energy transfer in the char zone. These reactions are valid up 

to 3000°F (1925°K). Above this temperature, additional reactions 
must be considered to accurately describe the reacting flow system. 

2. Under conditions of high mass fluxes ( > 0.05 Kg/irr-sec) 
and/or low temperatures (< 1365°K), the flow of pyrolysis products 
is essentially frozen. These conditions define the upper limitation 
of the frozen flow model. 

3. The. equilibrium flow model erroneously predicts the behavior 
in the char for all conditions except those in which the mass flux 
rate is smaller than 0.0001 lb/ft2-sec (0.0005 Kg/m2-sec). Mass 

flux rates greater than this value require the use of a non-equilibrium 
flow analysis to describe reactions occurring at finite rates. 


Reacting Flow of Pyrolysis Products Through Porous Graphite 

1. The same results are obtained for the reacting flow of pyrolysis 
products in porous graphite. The same important reactions and kinetic 
data that applied to the low density nylon-phenolic resin chars, 
likewise apply to the graphite between 1350 - 2300°F (1000°K - 1535°K) 
and 0.00003 - 0.10 lb/ft^-sec (0.00015 - 0.5 Kg/m^-sec). This 

permits the substitution of graphite for the brittle chars in 
studies requiring post- experimental analyses of the porous 
media; i.e., carbon deposition studies. 

2. In order to compare the results from char and graphite 
experiments, the mass flux in the pores must be identical. 

This corrects for differences in the porosity of each material. 

3. The inclusion of graphite "fines" in the pore spaces, 
during the fabrication procedure, results in reactions with 
hydrogen forming methane. Purging the graphite with helium 

or nitrogen prior to an experiment is sufficient to remove these 
dust- like particles from the pore spaces and eliminate the 
methane formed from this more reactive form of graphite. 
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Carbon Deposition Studies 


1. The use of radioactive tracers is a very effective 
method of determining the extent of carbon deposition within 
the char from the cracking of the pyrolysis gases. A plot of 
the percent radioactivity as a function of char depth very 
graphically determines the position where deposition starts , 
ends and is maximum. 

2. Carbon is deposited as a result of the thermal degra- 
dation of methane and phenol. Deposition starts near the middle 
of the char (z/L = 0.5) where the temperature is 1400°F (1035°K) 
and increases uniformly to a maximum value near the front surface 
at a temperature of 1900°F (1310°K). 


Degradation Product Distribution 

1. The products of methane degradation were determined 
by comparing radioactive tracer concentration profiles with the 
gas chromatograms of the exit gases leaving the Char Zone 
Thermal Environment Simulator. The species identified were 
carbon monoxide , carbon dioxide, ethylene, acetylene and 
unreacted methane. 

2. Similarly, the products of phenol degradation were 
determined. The species were identified as methane, carbon 
monoxide, carbon dioxide and unreacted phenol. 

3. The above methods provided excellent supporting 
evidence that the reactions considered important in the mathe- 
matical model were correct. 


Catalysts Evaluation Studies 

1. The effect of adding homogeneous catalysts (such as 
bromine) to the pyrolysis product stream to accelerate chemical 
reactions and, therefore, the heat absorption in the system 
was determined. Bromine (and the organic halides and other 
halogens) produced a catalytic effect which reduced the tempera- 
ture at which carbon deposition starts from 1400°F (1035°K) 

to 120Q°F (911°K). The position of maximum deposition was likewise 
shifted 300°F (150°C). A sharp decline in the carbon deposition 
probably indicated rapid depletion of carbon by reactions with 
water and/or carbon dioxide. 

The product concentrations of the exit gas stream is 
likewise different from the non-catalytic experiemnts and 
not predicted by the non-catalytic kinetic data. 

2. Molybdenum and tungsten heterogeneous catalysts did 
not affect a measurable change in the behavior of the system. 

A slight shift in the carbon deposition profiles were observed; 
however, the results were outside the experimental accuracy of 
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the analytical equipment, Platinum was used in earlier studies 
and, likewise, failid to produce a ©hang®. 


Alt Oxidation Studies 

The effect ©£ injecting air into low density nylon- phenolic 
resin chars at the front surface Is rapid erresien near the 
edge regions of the chars. These experiments were conducted by 
passing air into the char from the front surface with no 
counter flew of pyrolysis products from the rear. The 
experimental results were bracketed by the Isothermal analysis 
and values calculated using a linear description el the tempera- 
ture profile within the char. 


- 90 - 




NOMENCLATURE 


Symbol 

Description 

Units 

A 

area 

L 2 


species identification in a general expres- 
sion for a reaction 

none 


Runge-Kutta parameter in the energy equa- 
tion numerical solution method 

T 


ratio of system parameter defined in 
equation (2-13) 

none 

a 

formula number for chemical compounds 

none 


coefficients in the polynomial 
defined in equation (2-28) 

none 

B 

Runge-Kutta parameter in the species 
continuity numerical solution method 

M/L 2 t 


ratio of system parameter defined in 
equation (2-14) 

none 

b 

mass of a chemical element 

M 


coefficients in the polynomial defined 
in equation (2-29) 

none 

C 

ratio of the fluid to matrix thermal 
conductivity 

none 


molar concentration defined in 
equation (2-38) 

mole/L^ 

C 

P 

heat capacity of a pure component at 
constant pressure 

L 2 /t 2 T 

C 

P 

heat capacity of a mixture at constant 
pressure 

L 2 /t 2 T 
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Symbol, 

Description 

Units 

c 

total number of atoms of a chemical 
element In the system defined by 
equation (3-41) 

none 


coefficients In the polynomial defined 
in equation (2-35) 

non® 

D 

substantial derivative 

none 


mass diffusivity 

L 2 /e 

E 

energy of activation 

ML 2 /t 2 ' 

mole 

F 

free energy 

ML 2 /t 2 

F° 

free energy at a standard state of 298°K 
and 1 atom 

ML 2 /t 2 

f 

free energy function defined in 
equations (3-31) and (3-33) 

ML 2 /t 2 


any mathematical function 

none 

G 

augmented function of the quadratic 
approximation to the free energy 
function defined in equation (3-37) 

ML 2 /t 2 

g 

force of gravity 

L/t 2 

H 

enthalpy 

ML 2 /t 2 


enthalpy at a standard state of 
298°K and 1 atom 

ML 2 /t 2 

h 

heat transfer coefficient in equation 
(2-3) 

M/ t 3 t 


numerical integration increment sizes 

none 
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Symbol 

Descriotion 

Units 

J 

molar flux 

2 

moles/tL 

J 

mass flux 

M/tL 2 

K 

total number of gas species in the system 

none 

k 

thermal conductivity 

ML/t 3 T 

k° 

reaction rate constant 
frequency factor 

^function 
of react lot 
order. 

L 

char thickhess 

L 

M 

mobility of a porous matrix defined in 
equation (2-84) 

none 

M 

. w 

molecular weight 

M/moles 

ID 

total number of chemical reactions in 
the system 

none 

N 

molar flux of species in the char zone 

2 

mole/tL 

n 

moles of species 

moles 

n 

mole fraction of species 

none 

P 

pressure 

M/Lt 2 

P 

stoichiometric coefficient of the 
products in a chemical reaction 

none 

p' 

power on the concentration terms for 
products in a chemical reaction 

none 

Q 

volumetric flow rate 

1? It 


quadratic approximation of the free 
energy function (3-36) 

ML 2 /t 2 

q"' 

heat generation by chemical reaction 

M/t 3 
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Symbol 

Description 

Units 

q 

number of gas plus solid (condensed) 
phases in the system 

none 


energy transfer by conduction, 
convection or radiation 

M/t 3 

R 

ideal gas constant 

ML 2 /t 2 T 

mole 


chemical reaction rate 

* 

mole/tL* 

R 

effective chemical reaction rate for gas 
and solid species defined in equation (3-21) 

mole/tL" 

R 

e 

Reynold's number 

none 

R 

0 

dimensionless mole flux defined in 
equation (2-43) 

none 

r 

chemical reaction rate of a chemical species 

mole/tL’ 


stoichiometric coefficient of the 
reactants in a chemical reaction 

none 


radius of a cylindrical graphite test 
section included in equation (2-83) 

L 

r' 

power on the concentration terms for 
reactants in a chemical reaction 

none 

s 

power on the temperature in the rate 
equation (3-52) 

none 

T 

temperature 

T 

T' 

temperature gradient 

T/L 

t 

time 

t 

u 

ratio of the actual to estimated mole 
fraction in the free energy minimization 
calculation (3-42) 

none 

V 

volume 

L 3 
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Svmbol 

Description 

Units 

V 

velocity 

L/t 

W 

mass flux based on the superficial 
velocity 

M/L 2 t 

w 

p 

mass flux based on the velocity in the 
pores spaces of chars and graphite 

M/L 2 t 

X 

mass fraction of the species 

none 


distance in, the axial direction 
(Chapter II) 

L 

y 

distance in the radial direction (2-82) 

L 


estimated moles of species in the free 
energy minimization calculation 

none 

y 

estimated mole fraction in the free energy 
minimization calculation 

none 

z 

distance in the axial direction 
(Chapter II and III) 

L 

z 

compressibility factor 

none 

Greek 



“o 

ratio of the mass flux to the mass 
diffusivity used in equation ( 2-44 ) 

M/I. 4 

a 

viseous coefficient in the modified 
Darcy's "law equation (2-82) 

L- 2 

3 

inertial coefficient in the modified 
Darcy's law equation (2-82) 

L' 1 

Y 

permeability of a porous htedium 

l 2 

6 

Kronecher Delta 

none 

A 

a difference between two parameters 

none 

V 

del operator 

none 
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Symbol 

Description 

Units 

€ 

porosity of a porous medium 

none 

£ 

emissivity of a porous medium 

none 

i 

dimensionless char distance defined in 
equation (2-15) 

none 

8 

dimensionless temperatures defined in 
equations (2-11) and (2-12) 

none 

X 

parameter of a straight line (0-1) 

none 

n 

viscosity 

M/Lt 

ft 

3.1416 

none 

p 

density 

M/L 3 

a 

Stefan-Boltzmann constant 

M/t 3 T 4 

a 

collision diameter 

L 

T 

shear stress 

M/t 2 L 

0 

dimensionless parameters defined 
in equations (2-25), (2-26) and (2-27) 

none 


parameter used to calculate the viscosity 
of a gas mixture in equation (3-57) 

none 

£5 

collision integral 

none 

Ii 

Lagrangian multipliers in equation (3-37) 

none 

Subscripts 



c 

convection or conduction 


d 

carbonaceous deposit 


e 

effective or overall value 


f 

fluid or gas phase 
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Description 


Symbol 

g gas phase 

L front surface of a char 

m matrix or solid phase 

0 initial or inlet condition 

P pressure 

pyrolysis 

pores 

r radiation 

s solid phase 

T temperature 

total 

y distance in the radial direction 

z distance measured from the char back surface 

®° free stream condition 

1 inside surface of a cylindrical tube 

2 outside surface of a cylindrical tube 
Superscripts 

* a standard or reference condition 

° a standard or reference state 

1 derivative 
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